ASYMPTOTICS OF FORWARD IMPLIED VOLATILITY 



ANTOINE JACQUIER AND PATRICK ROOME 



(N 

O 
(N 



Abstract. We prove here a general closed-form expansion formula for forward-start options and the forward 
, implied volatility smile in a large class of models, including Heston and time-changed exponential Levy models. 

This expansion applies to both small and large maturities and is based solely on the knowledge of the forward 
characteristic function of the underlying process. The method is based on sharp largo deviations techniques, 
and allows us to recover (in particular) many results for the spot implied volatility smile. In passing we show (i) 
that the small-maturity exploding behaviour of forward smiles depends on whether the quadratic variation of 
the underlying is bounded or not, and (ii) that the forward-start date also has to be rescalcd in order to obtain 
non-trivial small-maturity asymptotics. 
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Ph ' 1- Introduction 



. 1-^ ^ Consider an asset price process (6''^')j>q with Xq = 0, paying no dividend, defined on a complete filtered 

probability space {^,J-, {J't)t>o,^) with a given risk-neutral measure P, and assume that interest rates are zero. 
In the Black-Scholes-Merton (BSM) model, the dynamics of the logarithm of the asset price are given by 



^ : (1.1) dXt = -l<J^dt + adWt, 

] where cr > is the instantaneous volatility and is a standard Brownian motion. The no-arbitrage price of 

■ the call option at time zero is then given by the famous BSM formula [T^llll]: Cbs('''7 k, cr) := E (c"^^ — e*^)^ = 
[ Af {d+) — e^-'A/" (d- ) , with d± := —-^j^ ± ^c-v/t", where M is the standard normal distribution function. For a 

given market price C°^^(t^W) of the option at strike e'' and maturity r we define the spot implied volatility 
(Tr (fc) as the unique solution to the equation C°'^*'(r, fc) ~ Cbs(t, fc, crr(fc))- 

For any t, r > and fc G K, we define [TOl |43] a Type-I forward-start option with forward-start date i, 

■ maturity r and strike e'^ as a European option with payoff (e'^'+^/e'^* — e'^)'*'. In the BSM model (jl.ip its 
^ ■ value is simply worth Cbs(''', fc, c). For a given market price C°'^^(i,T, fc) of the option at strike e*^, forward- 



start date i and maturity t we define the forward implied volatility smile oi^^ilti) as the unique solution to 
C°^^{t,T,k) = Cbs('''j ^, o't,r (fc)) since c^^CbsIt, fc, ct) > [TUl [M]- A second type of forward-start option 
exists [33] and corresponds to a European option with payoff (c''^'+^ — c'^+-^*)^. In the BSM model (|l.ip the 
value of the Type-II forward-start option is worth Cbs(''', k, a) [47]. Again, for a given market price C°^^'^^{t, t, k) 
of such an option, we define the Typc-II forward implied volatility smile at^rik) as the unique solution to 
C°^^'^^{t, t, k) = C^s{t, k, at.rik)). Both definitions of the forward smile are generalisations of the spot implied 
volatility smile since they reduce to the spot smile when i = 0. 

The literature on implied volatility asymptotics is extensive and has been studied using a diverse range of 
mathematical techniques. In particular, small-maturity asymptotics have historically received wide attention 
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due to earlier results from the eighties on expansions of the heat kernel [7]. PDE methods for continuous- 
time diffusions [5] , large deviations [TH [TB] , saddlepoint methods [H] and differential geometry [32] are among 
the main methods used to tackle the small-maturity case. Extreme strike asymptotics arose with the seminal 
paper by Roger Lee [12] and have been further extended by Benaim and Friz [5] [5] and in [301 [SB IMl US] • 
Comparatively, large-maturity asymptotics have only been studied in [51] HO] |36l |35j [22] using large deviations 
and saddlepoint methods. Fouque et al. |23| have also successfully introduced perturbation techniques in order to 
study slow and fast mean-reverting stochastic volatility models. Models with jumps (including Levy processes), 
studied in the above references for large maturities and extreme strikes, 'explode' in small time, in a precise 
sense investigated in [ll|2l|50l Ei |45l|18] . 

A collection of implied volatility smiles over a time horizon (0, T] is also known to be equivalent to the marginal 
distributions of the asset price process over (0,T]. Implied volatility asymptotics has therefore provided a set 
of tools to analytically understand the marginal distributions of a model and their relationships to market 
observable quantities such as implied volatility smiles. However many models can calibrate to implied volatility 
smiles (static information) with the same degree of precision and produce radically different prices and risk 
sensitivities for exotic securities. This can usually be traced back to a complex and often non-transparent 
dependence on model transitional probabilities or equivalently on the model generated dynamics of implied 
volatility smiles. The model dynamics of implied volatility smiles is therefore a key model risk associated with 
these products and any model used for pricing and risk management should produce realistic dynamics that 
are in line with trader expectations and historical dynamics. One metric that can be used to understand the 
dynamics of implied volatility smiles ([TU] calls it a 'global measure' of the dynamics of implied volatilities) is to 
use the forward smile defined above. The forward smile is also a market defined quantity and naturally extends 
the notion of the spot implied volatility smile. Forward-start options also serve as natural hedging instruments 
for several exotic securities (such as Cliquets, Ratchets and Napoleons; see [IHJ Chapter 10]) and it is therefore 
important for a model to be able to calibrate to liquid forward smiles. Despite the significant research on implied 
volatility asymptotics, there are virtually no results on the asymptotics of the forward smile: Glasserman and 
Wu |28| introduced different notions of forward volatilities to assess their predictive values in determining future 
option prices and future implied volatility, Keller- Ressel [40] studies a very specific type of asymptotic (when the 
forward-start date becomes large) , and empirical results have been carried out by practitioners in |10l I26[ I13j . 

We consider below a continuous-time stochastic process (Zj) and prove an expansion of option prices on 
(Zg) as e tends to zero. Setting = or = e^i/e then yields small or large- maturity expansions of 
option prices. This main result is presented in Section [5] as well as corollaries applying it to forward-start option 
asymptotics. We also translate these results into closed-form asymptotic expansions for the forward implied 
volatility smile (Type I and Type II). In Section [31 we provide explicit examples for the Heston, multi-Heston, 
Schobel-Zhu and time-changed exponential Levy processes. Section [4] provides numerical evidence supporting 
the practical relevance of these results and we leave the proofs of the main results to Section (5] 

Notations: Af{^J,,a'^) shall represent the Gaussian distribution with mean fi and variance cr^. Furthermore 
E shall always denote expectation under a risk-neutral measure P given a priori. We shall refer to the standard 
(as opposed to the forward) implied volatility as the spot smile and denote it CTt- The (Type-I) forward implied 
volatility will be denoted at^r as above. In the remaining of this paper e shall always denote a strictly positive 
small quantity. 
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2. General Results 

This section gathers the main notations of the paper as weh as the general results. The main result is 
Theorem 12.61 which provides an asymptotic expansion — up to virtually any arbitrary order — of option prices 
on a given process (X^), as e tends to zero. This general formulation allows us, by a suitable scaling, to obtain 
both small-time and large-time expansions. Indeed, setting = or = eXi/^ yields two expansions for 
different regimes, small and large maturities. The first rescaling is detailed in Section [2.1.11 and the second one 
in Section r2. 1.21 In each case, we shall make the computations explicit for the BSM p.ip case, which will also 
be needed to translate these expansions into expansions for the forward implied volatility in Section 12.21 Such 
expansions for European option prices and their corresponding spot implied volatilities are known for many 
models, and we shall consider here forward-start options (which clearly reduce to standard vanilla options when 
the forward-start date is null). 

2.1. Forward-start option asymptotics. Let {X^) be a stochastic process with re-normalised moment gen- 
erating function (mgf) 

fuX, 



(2.1) Ae(u):=£logE 



exp 



for all u G 2?^ 



where := {u E M. : |Ae (u) | < oo}. We now introduce the following critical assumptions. 
Assumption 2.1. For each u £ Vq q the re-normalised mgf can be represented as 

2 

(2.2) Ae (m) = A,,o (w) e' + O (e^) , as e tends to 0. 

i=0 

Further we suppose that for all e > the map Ag : 2?^ n- M is infinitely differentiable on Vq g C I?^ and e Vq g 
where we define I?g,o := {u G M : |Ag,o (u) \ < oo} and Vq q is the interior of I?o,o in K. 

The infinite differentiability assumption of the map A^ could also be relaxed by a C'*(I?g g) condition but 
this hardly makes any difference in practice and does, however, render some formulations awkward. If the 
expansion p.2|) holds up to some higher order n > 3, one can in principle show that both forward-start option 
prices and the forward implied volatility expansions below hold to order n as well. However expressions for the 
coefRcicnts of higher order are extremely cumbersome and scarcely useful in practice. 

Definition 2.2. [151 Definition 2.3.5] A convex function h :M.D Vh — > (—00,00] is essentially smooth if 

(i) is non-empty; 

(ii) h is differentiable in P^; 

(iii) h is steep, in other words lim„_j.oo |^'('Un)| = 00 for every sequence (Mri)neN in that converges to a 
boundary point of 2?,°. 

Assumption 2.3. Ag g is strictly convex and essentially smooth on g. 

Define the function A* : R R+ as the Fenchel-Legendre transform of Ag.g: 

(2.3) A*(fc) := sup {uk - Ag^o(ii)} , for all fc e R. 

For ease of exposition in the paper we will use the notation 

(2.4) A,,z(u) := aiA,,g(w) for / > 1, 
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where A^^o is defined in (|2.2p for i = 0, 1,2. The folfowmg lemma gathers some numediate properties of the 
functions A* and A,,- which will be needed later. 



Lemma 2.4. Under Assumptions lKl] and lKSl the following properties hold: 
(i) For any fc e R, there exists a unique u*{k) G q such that 

(2.5) Ao,i(ii*(fc)) = fc, 

(2.6) A* (fc) = u* {k)k - Ao,o {u* (fc)) ; 
(ii) A* is strictly convex and differ entiahle on R; 

(Hi) if a Cz Q such that Aoflla) = 0, then A*{k) > ak for all k G^\{AQ^i(a)} and A*(Ao,i(a)) = aAo.i(i)- 
Proof. 

(i) By Assumption l2.3l and l2.1l An i is a strictly increasing difFerentiable function from — oo to oo on I?o,o- 

(ii) By (i), 9fcA*(fc) — Aq i[(fc) for all fc S R. In particular 9fcA* is strictly increasing on R. 

(iii) Since Aq.i is strictly increasing, Ao,i(a) = fc if and only if u*{k) — a and then A*(Ao.i(a)) — aAo.i(a) 
using (|2.6p . Using the definition (|2.3p with a E T>q q and Ao,o(a) = gives A*(fc) > ak. Since A* is strictly 
convex from (ii) it follows that A*(fc) > ak for all fc G R \ {Ao,i(a)}. 



□ 



For ease of notation we shall write A^^; in place of Ajj (it*(fc)). Let / : R_|- R_|- be a function such that 
(2.7) /(e)^ = c + 0(e), for some c > 0, as e tends to zero. 

For any 6 > we now define the functions Af,, : R \ {Ao.i(O), Ao.i(&)} x (0, oo) R by 



(2.8) .n^..._bV^k^>^,+ey^fie)l,,.o, 

u*{k){u*{k)-h).JWA^2 ' 

(2.9) ^(fc,s) :=c"^-W/^+'=/(^'A(fc,s) (^l + T{b,k)e+'^^^0^1^^^1{,^,} 



u*{k) 



where T : [0, oo) x R\{Ao,i(0), Aoa(6)} 



(2.10) 



T(&,fc-) := A2,o- 



24A3.2 



i is given by 
4Ai,iAo,3 + An,4 A? , + Ai ■ 



Ao,: 



46=0} 



Ao,3 



2A, 



0,2 



2u*(fc)A2 2 2(it*(fc)-6)A2, 



Ai,i (6- 2u*(fc)) +3 



U*{k) {U*{k) - b) Ao,2 U*{kf {U*{k) - by Ao,2 ' 

Remark 2.5. The domain of definition of excludes the set {Ao4(0), Ao4(6)} = {fc G R : M*(fc) € {0,6}}. 
For all k in this domain, Ao,2('«*(fc)) > by Assumption 12.31 so that Ab is a well-defined real- valued function. 

The main result of the section is the following theorem on asymptotics of option prices. A quick glimpse at 
the proof of Theorem 12.61 in Section [5TT] shows that this result can be extended to any arbitrary order. 
Theorem 2.6. Let {X^) satisfy Assumvtions \2.1\ and \2.3l and f : R+ R+ be a function satisfying (|2.7p 
with constant c G Vq q H R-|-. Then the following expansion holds for all k G R\{Ao,i(0), Ao,i(c)} as e \ 0; 

+" 

l{fe>Ao,i(c)} 



Ac{k,e) = E 

-E 



Hfe<Ao,i(0)} 



l'-{Ao,i(0)<A:<Ao.i(c)}' 



where Aq.i is defined in (j2.4p and Ac{k,e) in (j2.9 
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Remark 2.7. In the case c = 0, the expansion holds for all k G R\{Ao,i(0)} and the last term on the right-hand 
side disappears, since the indicator function is taken over an empty set. 

Let {Xt)t>o be a stochastic process. For any t > 0, we define (pathwise) the process 
(2.11) Xt+r - Xt 

We now specialise Thcorem l2.6l to forward-start option asymptotics. 



2.1.1. Diagonal small-maturity asymptotics. We first consider asymptotics when both t and r are small, which 
we term diagonal small-maturity asymptotics. 



Corollary 2.8. // ( xif ^ 



e>0 



satisfies Assumvtions lKli \ 2.S\. then the following expansion holds: 



-A*(fe)/e+fc+Ai,0g3/2 



T(0,fc) 



u*{k) 



■E 



(e'^-e 



X 



(et) 



4fe>Ao,i(0)} 



IL{fc<Ao.i(0)}, 



as £ tends to zero, where Ao,i, A*, u*{k), T and Aj,/ are defined in ^T^, (P3)) . (|2.10p and (HH]). 

Proof. Set {Xfr) := (xir*'' ) and / = 1. Then c = and the corollary follows from Theorem 



□ 



Corollary 2.9. In the BSM model (|l.ip the following expansion holds as e tends to zero: 



E 



l{fc>o}+E 



(e^-e^i^^'')" 



4fe<o} 



3/2 



3 1 

fc2 ^ 8 



Proof. For the rescaled (forward) process ( xi^*' ) in the BSM model (11. ip we have KJu) = Aq o(it)+£^i o(u) 

V /e>0 

for li e M. where Ao,o(w) = u^cr^r/2 and Ai^o(w) = —ua'^T/2. It follows that Ao,i(m) = ucr'^r, Ao.2('«) = cr^r and 
Ai,i(u) = — o'^t/2. For any A; G M, u*(fc) := k/{a'^T) is the miique solution to the equation Ao,i(u*(fc)) = k and 



A*(fc) = k /(2a t). Aq.o is essentially smooth and convex on R and the BSM model satisfies Assumptions 12.1 
and 12.31 Since e Pqq and Aq.i(0) = 0, the result follows from Corollary 



□ 

It is natural to wonder why we considered diagonal small-maturity asymptotics and not the small-maturity 
asymptotic of at^r for fixed t > 0. In this case it turns out that in many cases of interest (stochastic volatility 
models, time-changed exponential Levy models), the forward smile blows up to infinity (except at-the- money) 
as T tends to zero. However under the assumptions given above, this degenerate behaviour does not occur in 
the diagonal small-maturity regime (Corollary 12. 8p . We leave the precise study of this degeneracy for future 
research, but provide a preliminary conjecture explaining the origin of this exploding behaviour. Consider a 
two-state Markov-chain dXt = —^Vdt + \/V dVFt, starting at Xq = 0, where is a standard Brownian motion 
and where V is independent of W and takes value V\ with probability p g (0,1) and value V2 <E (0,Vi) with 
probability \ ~ p. Conditioning on V and by the independence assumption, wc have 



E 



for all u E 



Consider now the small-maturity regime where e = r, f{e) — 1 and X^ xi'^ for a fixed t > 0. In this case 
an expansion for the re-scaled mgf in p.2|) as r tends to zero is given by 



Ae(u) = TlogE (^c"(^*+--^"*)/^) = + rlog (pc-^^^/^) + tO 



tO c 



-u^(Vi-V2)/i2 



for all u e 
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Since Vi > Yi the remainder tends to zero exponentially fast as t tends to zero. The assumptions of Theorem l2.6l 
are clearly satisfied and a simple calculation shows that limT-^^o <^t.T{k) = \fV\. This example naturally extends 
to an n-state Markov chain, and a natural conjecture is hence that the small-maturity forward smile does 
not blow up if and only if the quadratic variation of the process is bounded. In practice, most models are of 
unbounded quadratic variation (see the examples in Section [3] below), and hence the diagonal small- maturity 
asymptotic is a natural scaling. 



2.1.2. Large-maturity asymptotics. We now consider large-maturity asymptotics, when t is large and t is fixed. 
Define the function : R\{Ao,i(0), Ao4(l)} x (0, oo) ^ M by 



(2.12) 



r(A'(fc)-fc)+Ai. 0^.-1/2 



'{k) {u*(k) - 1) ^2^ 



where A*, u* , T, A,j are defined in (gS]), ([23]), (j^ . From Remark[231 the function B°° is well-defined. 

Wc now have the following large-maturity asymptotic for forward-start options. 

Corollary 2.10. If the process (r^^Xr^j satisfies Assumvtions \2. ll \2.S\ with e = r^^ and 1 G T^qq, then 



r>0 



the following expansion holds as r tends to infinity: 



E 



(e^^*' - e^- 



l{fc>Ao,i(l)}+IE 



4fc<Ao,i(0)}- 



e^^^'Ae^- 



l{Ao,i(0)<fc<Ao,i(l)} = B°°{k,T), 



where Aq.i is defined in (|2.4p and B°° in (|2.12p . 



Proof. Let {X^) := (^eATj^^j , e := l/r and /(e) = 1/e (c = 1), then the result follows from Proposition 12. 61 □ 
In the BSM case dTT]), define the function B^g : M\{-ctV2, a'^ /2} x (0, oo) ^ M by 



(2.13) i?g^s(fc,T):=exp -r 



(fc + aV2)' 
2a2 



k 



[Ak^ - cr4) ^ 



4ct2 -t- 12A:2; 
(4P-(t4)V 



and we have the following corollary. 

Corollary 2.11. In the BSM model (jl.ip the following expansion holds as r tends to infinity: 



E 



(e^^"-e'^- 



l{fe>o-V2} 



l{-crV2<fc<<TV2} + 



kT Y**) 



Proof. Consider the process ix'^^ It] and set e ~ t ^. In the BSM model (11. ip for any w G M. we have 

V /t>0 



Ae(u) = logE(exp(wX| 0) = Ao,o(w) = icr2y(u - 1). Thus Ao,i(u) = a^{u - 1/2) and Ao,2 (w) = (t^ For 
any fc G R, the equation Ao^i(w*(fc)) = fc has a unique solution u*(fc) = 1/2 + k/a"^ and therefore A*(fc) = 
[k + 0-2/2)^ /(2cr2). Aq.o is essentially smooth and strictly convex on R and Assumptions 12 . 1 1 and that 12.31 are 
satisfied. Since {0, 1} C I?o o ^li*^ result follows from Corollarv l2.10l □ 



2.2. Forward smile asymptotics. In this section we translate our results on forward-start options into asymp- 
totics of the forward implied volatility smile k ^ at.T (k). We first focus on the diagonal small-maturity case. 
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For i = 0, 1, 2 we define the functions : M* x [0, oo) x (0, oo) — > R by 

fc2 



voik,t,T) 



2TK*{k)'- 



2TUo(fc,<,T)2 / ^2gAi,o(n*(fc)) \ 

vi{k,t,T) := — log' 



(2.14) 



U*ik)^^Ao^2iu*{k)){TVo{k,t,T)f\ 



H , . - j^vo{k,t,T)vi(k,t,T), 

vo[k,t,T) k^ 

wliere A*, w*, A^,;, T are defined in (gS]) ([211), (gUHl)- Also if Ao.i(O) = then A*(fc) > for k e R* and 

A*(0) = from Assumption 12.11 and Lemma [2.4^ 11) so that vq is always strictly positive. A direct application 
of L'Hopitals rule, together with Lemma r2.4f i')fii). Assumption 12.11 and Ao,i(0) = 0, shows that for any t > 0, 
T > 0, the map vo{-,t,T) can be extended by continuity at the origin with vo{0,t,T) = 1/(tu* (0)). All the Vi 
{i — 0,1,2) are hence well-defined real-valued functions (see also Remark 12. 5p . The diagonal small-maturity 
forward smile asymptotic is now given in the following proposition, proved in Section 15.11 

Proposition 2.12. Suppose that [Xe%^^\ satisfies Assumvtions \2. 1\ and \2.'A and that Ao i(0) = (defined 
in (|2.4p ). The following expansion holds for the corresponding forward smile for all fc G M* as e tends to zero: 

(2.15) cret,erik) = wo(fc, t, t) + vi{k, t, T)e + V2ik, t, T)e^ + O (e^) , 

where vq , vi and V2 are given in p.l4p . 

In the large-maturity case, define for i = 0, 1, 2, the functions : M\{Ao.i(0), Ao.i(l)} x [0, oo) — > M by 



v'^{k,t) 



2 2A*(fc) -k- 2^K*{k){K*{k) - k) , if fceM\[Aoa(0),Ao,i(l)] 



2 2A*(fc) - fc + 2v/A*(fc)(A*(fc) - k) , if fcG (Aoa(0),Ao,i(l)), 



oon .^ 8w§°(/c,t)2 / / 4fc2_„g°(fc,t)2 

^T{k,t) :=^U-L^ Ai,o("*(fc))+log' 



(2.16) 



4k-^-v^[k,tY \ \4{u*{k) - l)u*{k)v^{k,tf/^^Ao,2{u*{k)) J J ' 

vf(k,t) := ^ \8k'^v?° (fc, t)v^ {k,tf ivf (fc, t) + 6) - 16fcS?° ik, tf 

v^{k,t){v^{k,tf-ik^f^ 1 : ; V , ^ V 1 V , ; 

- 2T(1, k)v^{k, tf (vg^{k, tf - 4k^^ ' - k^v^ {k, tf (§6 + v'{°{k, tf + Sv'^ (fc, i)) 

-v^{k,tf{vr{k,t) + 8) . 

A* is defined in [131 u* in dSS]), A^,, in (EH) and T in (PJU)) . Since {0, 1} C V^^ and Ao,o(l) = Ao,o(0) = 0, 
we always have A*{k) > max(0, fc) from Lemma [2.4( 111'). One can also check that < v^{k,t) < 2\k\ for 
k e M\ [Ao,i(0), Ao,i(l)] and v^{k,t) > 2\k\ for k e (Ao,i(0), Ao,i(l)). Together with Remark[231 this implies 
that the functions (i = 0, 1, 2) are always well-defined and real- valued. 

Remark 2.13. By Assumption 12.11 and Lemma [2.4( 111) we have A*(Ao,i(0)) ~ 0. Further by the assumptions 
in Proposition 12.141 below we have 1 G Vq q and Ao.o(l) = 0. Again from Lemma [2.4f iii) this implies that 
A*(Ao,i(l)) = Ao,i(l). Hence for all t>0, v^{-,t) can be extended by continuity on M. 

The large- maturity forward smile asymptotic is given in the following proposition, proved in Section 15.11 
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Proposition 2.14. Suppose that ( t^^xI*"* J satisfies Assumvtions \2. l\\2.S\ with e = and that Aq q{1) = 

Ar>0 

with 1 e g (all defined in A ssumption \2.1\) . The following expansion then holds for the forward smile as r 
tends to infinity for k € K\{Ao,i(0), Ao,i(l)}.' 

(2.17, + ^ + ^ + 

where , and are defined in (j2.16p . 
Remark 2.15. 

(i) If we set t = in (|2.15p and (|2.17p then we recover — and actually improve — the implied volatility asymp- 
totics obtained in [H], [ig, ES], [H], [22. 

(ii) The forward smile results can be extended to a deterministic interest rate setting by considering the 
forward price instead of the stock price and re-scaling the strike appropriately. 

For the (J"tj)-martingale price (e'^")„>o (under P) define the stopped process X* := Xt^u for any i > 0. 
Following [43] define a new measure P by 

(2.18) F{A) := E {e^'^+^l^ , for every A G Ft+r- 

The stopped process (c^")„>o is a (J^tAii)ti-martingalc and ()2.18p defines the stopped-share-price measure P. 
The following proposition shows how the Type-II forward smile at.r can be incorporated into our framework. 

Proposition 2.16. // (e^'')(>Q is a (Tt) -martingale under P, then Propositions \2.12\ and \2.l4\ hold for the 
Type-II forward smile at^r with the mgf (|2.ip calculated under f. 



Proof. We can write the value of our Type-II forward-start call option as 



E 



E 



E 



E 



Proposition l2.6l and Corollaries l2.8[[2.10l hold in this case with all expectations (and the mgf in (|2.ip ) calculated 
under the stopped measure P. An easy calculation shows that under P, the forward BSM mgf remains the same 
as under P. Thus all the previous results carry over and the proposition follows. □ 



3. Applications 

3.1. Heston. In this section, we apply our general results to the Heston model, in which the (log) stock price 
process is the unique strong solution to the following SDKs: 

dXt =-lvtdt+^tdWt, Xo = 0, 

(3-1) dVt K{e -Vt)dt + CVVtdZt, Vo^v>0, 

d{W,Z)^ ^pdt, 

with K>0, ^>0;6'>0 and \p\ < 1. The Feller SDK for the variance process has a unique strong solution by 
the Yamada-Watanabe conditions [351 Proposition 2.13, page 291]). The X process is a stochastic integral of the 
V process and is therefore well-defined. The Feller condition, 2k.O > ensures that the origin is unattainable. 
Otherwise the origin is regular (hence attainable) and strongly reflecting (see [351 Chapter 15]). We do not 
require the Feller condition in our analysis since we work with the forward mgf of X which is always well-defined. 
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3.1.1. Diagonal Small- Maturity Heston Forward Smile. The objective of this section is to apply ProDOsition l2.12l 
to the Heston forward smile. We define the function A : Vt ^- x [0, oo) x (0, oo) — M by 

uv 

(3.2) A (u, t, t) ■= — — . _ r r rPn~' ^'^^ " ^ 

where 

(3.3) I?t,r := |u e M : A(u,0,t) < ^1 and p := \/l~p^. 
Further we let the function L : Vt.r x [0, oo) x (0, oo) —s- M be defined as 

(3.4) L{u,t,T) := Ln{u,T) + k{u,t,Tf [ /a - ^ — - KtA(u, t, r) - — ^ log 1 - ^ , ; 



where the functions Li : T>t^T x (0, oo) — > R for i = 0, 1, are defined by 
(3.5) 

io(it,T):=7j (it/3-ao)i™-21og' 

\ \ ^ .yu / / 



i±(p - do) (1 - e-^'*"-') (.91 - idigoru) 



(1 - ,00-^^0-) (^^^^ " + - - ^ ' l-,oe--.- 

with 

, , i(2Kp-?) ip-p , 2k -^p 

do := Cp, «i := izz , 90 ■= -. — — and gi := — 

2p ip + P ^p{p+ip) 

Remark 3.1. For any t > 0,t > the functions Lq and Li are well-defined real-valued functions for all u G I?t,r 
(see Remark [5?6l for technical details). Also since A(0, t)/A(0, 0, r) = 1, L is well-defined at m = 0. 

Proposition 3.2. In Heston, Corollarv \2.8\ and Proposition \2.12\ hold with I?o,o = 'Dt.r, Ao,o = A and Ai_o = L. 

Proof. We simply outline the proof of the proposition, and we refer the reader to Section [5. 2. II for the details. 

(i) In Lemma [5751 we show that 2?o.o = T^t,T and G q; 

(ii) In Lemma 15.51 we show that the Heston diagonal small-maturity process has an expansion of the form 
given in Assumption 12 . 1 1 with Kq q = A and Ai q — L; 

(iii) In Lemma 15.71 we show that A is strictly convex and essentially smooth on I?° ^ , i.e. Assumption 12.31 

(iv) Ae is infinitely differcntiable and Ao.i(O) = 0. 

We now apply Proposition 12.121 and this completes the proof. □ 

In order to gain some intuition on the role of the Heston parameters on the forward smile we expand (j2.15p 
around the ATM point in terms of the log strike k. We now define the following functions: 

,yo{t, r) := ^ (24^0 + e (p' - 4) + I2v{ip - 2n)) -\{e+ {v - 9)) , 
Mt, r) ^ {e (1 - p') -2n{v + 0)+ ^pv) + 



er 



(3.6) 

V2it, t) := (8Ok0 (I3p2 - 6) + (521^^ - 712^^ + 176) + AOp^v {^p - 2k) 

- ll^ (4-^ (16 - 7p') + {7p' - 4) {9e + 4«.) ) + ^ (v 
The proof of the following corollary is given in Section r5. 2. II 
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Corollary 3.3. The following expansion holds for the Heston forward smile as e and k tend to zero: 
(3.7) <j',,(k) = V + sMt, r)+(4+ r)) k + f ^^-^^ + |^ + eu^it, r)] k' 



+ 0{k^) + 0{ek^) + 0{e^). 

Remark 3.4. The following remarks should convey some practical intuition about the results above: 

(i) For t = this expansion perfectly lines up with Corollary 4.3 in [21] for the implied volatility smile. 

(ii) Corollary 13.31 implies cret,er(0) — cro,eT(0) = — (^^ + 4k(w — 0)) + 0{e'^), as e tends to zero. For small 
enough maturity, the spot ATM volatility is higher than the forward one if and only if + 4k(u — 9) > 0. 
In particular, when v > 9 the forward ATM volatility is lower than the corresponding spot ATM volatility 
and this difference is increasing in the forward-start dates and volatility of variance. In Figure [2] we plot 
this effect using 9 ^ v and 9 > v + (4k). The relative values of v and 9 impact the level of the forward 
smile vs spot smile. 

(iii) Similarly, we can deduce some information on the forward skew from CoroUarv 13.31 

f) ^ (0)- (4t/i(t,r)z;-epi^o(t,T)) ^ 



and hence 



, . ^pt (3£.^ + 4:k(v - 9)) , „ 



(iv) Likewise an expansion for the Heston forward convexity as e tends to zero is given by 

2 _ e^((2-5p^)r + 60 v^{t, T)g^(3t + (1 - ^p'Y) + 6r^(pC^i(f, r) - ^^^{t. t)v) ^ , 

and in particular d^a^t.ET (0) — 9^cro,er (0) = ^^t/(4r7j^/^) + 0{e). For fixed maturity the forward convexity 
is always greater than the spot implied volatility convexity (see Figure [2|) and this difference is increasing 
in the forward-start dates and volatility of variance. At zeroth order in e the wings of the forward 
smile increase to arbitrarily high levels with decreasing maturity, (see Figure [IJa)) This effect has been 
mentioned qualitatively by practitioners [T3|. As it turns out for fixed t > the Heston forward smile 
blows up to infinity (except ATM) as the maturity tends to zero. This is clearly outside the scope of our 
main theorem, and we leave this degenerate case for future research. 

In the Heston model {e^*)t>o is a true martingale [H Proposition 2.5]. Applying Proposition 12.161 with 
Lemma 15.81 giving the Heston forward mgf under the stopped-share-price measure, we derive the following 
asymptotic for the Typc-II Heston forward smile cTt,T- The proof of CoroUarv 13.51 is omitted as it is analogous 
to the proofs of Proposition [3?2] and CoroUarv 13.31 Set 



Mt,T) ■.= l^o{t,T)+^pvt, Vl{t,T) ■.= Mt,T), V2{t,T) -.^ ^^{t, t) + ^ {7 ^ i) - , 

4:8v ^ ' 8vT 

with vq, vi and V2 defined in l3.6l In particular when p = or t = 0, = t'i (i = 1, . . . , 3), and hence the Heston 
forward smiles Type-I and Type-II arc the same as shown in the following corollary. 

Corollary 3.5. The diagonal small-maturity expansion of the Heston Type-II forward smile as e and k tend to 
zero is the same as the one in Corollarv \3.3\ with vq, vi and replaced by Vq, vi and v^. 
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(a) Small-maturity forward smile explosion. 



(b) Type I vs Type II forward smile. 



Figure 1. In (a) we plot forward smiles with forward-start date t = 1/2 and maturities 
T = 1/6,1/12,1/16,1/32 given by circles, squares, diamonds and triangles respectively using 
the Heston parameters v ~ 0.07,6 ~ 0.07, k ~ l,p ~ —0.6,^ = 0.5 and the asymptotic in 
Proposition 13.21 In (b) we plot the Type I (circles) vs Type 2 (squares) forward smile with 
t = 1/2, r = 1/12 and the Heston parameters v = 0.07,6* = 0.07, k = l,p ^ -0.2, ^ = 0.34 
using Corollaries 13.31 and 13.51 




Figure 2. Forward smile vs spot smile with v = 9 and 6 > u + $^/ (4k). Circles (t = 0, r = 
1/12) and squares {t = 1/2, r = 1/12) use the Heston parameters v ^ 6 = 0.07, k = 1, p = —0.6, 
^ = 0.3. Diamonds {t = 0,t = 1/12) and triangles {t = 1/2, t = 1/12) use the same parameters 
but with 6* = 0.1. Plots use the asymptotic in Proposition [ 



3.1.2. Large-maturity Heston forward smile. We apply here Proposition 12.141 to the Heston forward smile. We 
shall use the standing assumption k > p£^, needed in the proof of Proposition 13.81 If this condition fails then 
we have a finite explosion time for moments greater than one for the price process (|3.ip and consequently the 
limiting mgf is not essentially smooth on its effective domain and Assumption 12.31 is violated. This a standard 
assumption in the large- maturity implied volatility asymptotics literature [20l [22l [36] . It has been relaxed in [35] 
for the Heston model to study the limiting spot smile, but not the higher-order terms. This restriction bears no 
consequences in markets where the implied volatility skew is downward sloping, such as equity markets, where 
the correlation is negative. Define the quantities 



(3 



u± 



and 



2e(e« 



1) 
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with 



(3.9) 77 := ^/e(l-p^-) + (2k-pC)2, „ ■= _ 16^26''* and V := ^(e"* - 1) - iKpe"*, 

as well as the interval I?oo C M and the real numbers p_ and p+ by 

7i!j_) , if — 1 < p < p_ and t > 0, 

, if /9+ < p < min(l,K/^),t > and k > p+£,, 
[u_ , u+] , a p- < p< min(p+ , k/^) , 

g-2«t |'^(g2Kt _ _^ ^ l)^16K2e2«t +^2(1 _(,,.t)2 

(3.11) p± := ^ 



8k 

Remark 3.6. The following remarks are proved in Lemmas 15. 9[ 15.101 and ISTT^ and we summarise them here. 

(i) When k/^ < the second case in p.lOp never occurs. From the proof of Lemma [5.10( 1) if k > ^ then 
p+ < 1 so that miii(/9_|_, k/^) < 1. 

(ii) If t > and p < p^, then u+ > u*^ > 1 and if i > and p > p+, then u_ < u*_ < 0. 

(iii) We always have —1 < p- < and if k > then 1/2 < p+ < 1. Also = —1 if and only if t = 0. 
When K > p^^ then = 1 if and only if t = 0. Finally 1/ defined in (j3.9p is a well-defined real number 
foraUpG [-l,p_]U[p+,l]. 

We define the functions V and H from Deo to R by 

(3.12) F(u) := — (k - p^M - d(u)) and H (u) :^ —^^-^——^ - — log ' ^* ^ 



Ke* - 2/3tT/(M) ^2 \^k6I(1-7(u)) 
with 

(3.13) d{u):={{^-pCur + ue{l-u)f^\ ^(^u) := and A := ^ (l - c^*) . 

Remark 3.7. Wc have (k - p^u)^ + u (1 - it) ^2 > q and k9 - 2l3tV{u) > for all u G I?oo from the proof of 
Proposition 15.131 Further by the definition of 7 in p.l3p we have 7(7/) E (—1, 1) for all u G 2?oo using (|5.38p 
and (|5.39p in the proof of Proposition 15. 131 So V and H are always well-defined real- valued functions. 

Finally we define the functions 5* : K — > [u_ , it+] and y * : R — > R+ by 

^2c2 , o^..fl^f ^ ,.2/12^-1/2 



(3.14) g*(x) := ^ '- and V* (x) := q* {x)x - V (q* (x)) . 



^ - 2kP + {K0p + X^) 77 (x2^2 _^ 2xK0p^ -f- k2 

2^(1 

where ^ and 77 are defined in p.l2p and in p. 91) . The following proposition gives the large- maturity forward 
Heston smile in Case (iii) in p.lOp . 

Proposition 3.8. Suppose that P- < p < min (p_|-, k/^) in the Heston model with p± defined in (j3.11l) . Then 
Corollary \2.10\ and Proposition \2.1^ hold with Aq.q — V , K* = V* , u* ~ q* , Ai g = H, A2_o = and 
T^o.o [^^-7""+]; where V, H, V* , q* and u± are defined in p.l2p . p.l4p and p.Sp . 

Proof. We simply outline the proof of the proposition, and we refer the reader to Section r5. 2. 21 for the details. 

(i) In Proposition 15. 131 we show that I?o,o = T^oo and that {0, 1} C I?^; 

(ii) In Lemma 15.141 and Remark 15. 161 we show that the process has an expansion of the form given in Assump- 
tion [2lT] with Ao.o = V, Ai.o = H and A2 = 0; 
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(iii) By Proposition 15.131 and Lemma 15. 9[ V is strictly convex and essentially smooth on if /9_ < p < 
min (p+j /t/^); see also Remark [XW ii) : 

(iv) Ae is infinitely differentiable and V{1) = from Lemma \UM 

(v) u* can be computed in closed- form and is given by q* in p.l4p . 

A direct application of Proposition 12.141 completes the proof. □ 

Remark 3.9. 

(i) In the Hcston model there is no t-depcndence for in (|2.17p since V* docs not depend on t. Therefore 
under the conditions of the proposition, the limiting (zerotli order) smile is exactly of SVI form (see [27]). 

(ii) For Cases (i) and (ii) in p.lOp the essential smoothness property in Assumption 12.31 is not satisfied and a 
different strategy needs to be employed to derive a sharp large deviation result for large-maturity forward- 
start options. We leave this analysis for future research. 

(iii) For t = we have p± = ±1 and Proposition l3. 81 agrees with and extends the Heston large-maturity implied 
volatility asymptotics in [50] and [22] • 

(iv) The condition p G [0, min(l/2, k/^)] is stronger than the condition in Proposition [2111 

(v) Even though the rate function V* does not depend on t, there is t-dependence through p± and the 
function H (see the ATM example below). That said, to zeroth order and correlation close to zero, 
the large-maturity forward smile is the same as the large-maturity spot smile. This is a very different 
result compared to the Heston small- maturity forward smile, as mentioned in Remark 13. 4f iv). where large 
differences emerge between the forward smile and the spot smile at zeroth order. 

We now give an example illustrating some of the differences between the Heston large-maturity forward smile 
and the large- maturity spot smile due to first order differences in the asymptotic (|2.17p . This ties in with 
Remark I3.9r v'). Specifically we look at the forward ATM volatility which, when using Proposition 13.81 with 
P- < p < min(p+,K/^), has the asymptotic 

0-2^(0) = v^{0) + v'^{0, t)/T + O (l/r^) , as t tends to infinity. 



with 



«o°°(0) 



AeK{i] - 2k + fp) 



e (1 - P^) ' 

^ _ l&Kv{pi-2K + 7^) I&kO ( Ae-«* (2k - + (l - 2^^) t^) 

^1 1^1, I) — ^ 1 iOg 



e 8K(i-p2)2^ 

Slog 



C(l-p2)'^V^(2ep-4K + 2r;) 



(ai - 2p2) - p(r? - 2k)) {p{^ - 2k) + ; ' 

77 is defined in p.9p and A := 2k (l -t- e''* (l - 2^^)) - (1 - e''*) (p^ + 77). To get an idea of the t-dependence of 
the ATM forward volatility we set p = (since Proposition 13 . 81 is valid for correlations near zero) and perform 
a Taylor expansion of f^(0, t) around t = 0: 

When V > 6 then at this order the large r-maturity forward ATM volatility is lower than the corresponding 
large r-maturity ATM implied volatility and this difference is increasing in t and in the ratio f /k. This is similar 
in spirit to Remark 13.4^ 11) for the small-maturity Heston forward smile. 
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3.2. Multivariate Heston. The n-Heston model (n e N) is defined as tlie unique strong solution to the 
following SDE: 

n n . 



2 

i=l 



(3.15) d^/'"' = AC,; (0, - f/'"' ) dt - ■^/''dzf ' , z = 1, . . . , n Fq^*' = 4'' G 



^+1 



d(W^W,H^O))^ =d(zW,Z«)^=l{,^,}dt, j,j = l,...,n, 

where /« = (ki,...,k„) € (0,oo)", p G (-1,1)", 6 G (0,oo)", | G (0,oo)", t^o G (0,cx))". The independence 
assumption of the variances allows the forward mgf of the n-Heston to be written as the sum of the forward 
mgf's of the individual Heston models. The asymptotics of the forward smile are then given directly in terms 
of the results in Section l3.ll For ease of notation in this section, whenever we reference a function or variable 
used in the Heston analysis in Section [3. II and use the index i, it means that function or variable defined using 
the i-th Heston model parameters in p.l5p . So for example when we reference the function Li in p.4p we mean 
that function L evaluated using the i-th Heston model parameters in p.l5p . 

Proposition 3.10. In the n-Heston model Corollary \2.8\ and Proposition \2.12\ hold with 2?o,o = nr=i -^tr) 
Ao,o = J2i=i ^''^'^ Ai.o = J2i=i where T>1.^, Ai and Li are defined in \3.3\ \3.2\ and \3.4\ 

Consider for instance n — 2, and define the functions: 

2 

(3.16) i^oit, t) ^ [( V 240,K, («i + v^f + I2v, («i + v^f {^^p, - 2k,) 

48 {V1+V2) 

+ {p1 - 4) [v, + V2f ) - vl^l {pI - 4) - vlil {p\ - 4) - 2v^V2 {H (5p? - 2) 



- 966P1P2 + il (5p2 - 2) j J - ^^^^^^^^^ 4AC, (^;i + v^) {v, - Oi) + j , 
(3.17) Mt, r) := -— ^ 3 f ^ (4 - 7p^) vj + 2v,V2 {^pl + 2) - 156^1^1^2 + ^| {^pl + 2)) 



t 



\i=l 
2 



4r (wi + 



,1=1 / 

In order to gain some intuition on the role of the Heston parameters on the forward smile we expand our solution 
around the ATM point in terms of the log strike k. 

Corollary 3.11. The following expansion holds for the 2-Heston forward smile as e and k tend to zero: 

al ,,{k) = +V2+ eMt, r) + ^lEl^ll^P^^k + vi {t, t)^ + ©(fc^) + 0{ek) + ©(e^). 

Remark 3.12. Remarks similar to Rcmark l3.4f ii')-fiv) for the Heston forward smile also apply to the 2-Heston. 

Proposition 3.13. If < pi < min(/9!|_, k/^) for i = 1,2,..., 71 (p^. defined in ()3.1ip ). then Corollary \2.10\ 
and Proposition \2.14\ hold with Ao,o = ^"=1 ^i-O = ^7=i ^2,0 = and I?o,o ~ 0"=: '^V] ; with Vi, 
Hi, defined in (|3.12p and (|3.8p . 
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3.3. Schobel-Zhu. The Schobel-Zhu (SZ) stochastic volatihty model [48] is an extension to non-zero correlation 
of the Stein & Stein [33] model in which the logarithmic spot price process (^t)t>o satisfies the following system 
of SDEs: 

1 9 

dXt = --cTt^t + (JtdWt, Xo = xoe M, 

(3-18) dat =n{0-at)dt+-^dZt, ao = 07 > 0, 

d{W,Z)^ =pdt, 

where n, 9 and ^ are strictly positive real numbers, p G (—1, 1) and (W^t)t>o (-^Ooo ^-""^ standard 
Brownian motions. The volatility process (o't)j>o is Gaussian and hence the SDE is well-defined. The process 
{Xt)fyQ is simply the integrated volatility process and hence is well-defined as well. The analysis in this section 
is very similar to the diagonal small-maturity Heston analysis and therefore the proofs are omitted. Note that 
although in some cases we use the same variables as in the Heston analysis they may have a different definition 
in this section. We present limited results to highlight the similarities and differences between the Heston and 
Schobel-Zhu forward smiles. 

Proposition 3.14. In the Schobel-Zhu model Corollary and Provosition hold with T^q q = 'Dt,T OLi^d 
Ag = A, where Vt ^ and A are defined in \3.3\ and[ 



Remark 3.15. At zeroth order in e the SZ diagonal small-maturity forward smile is the same as in Heston 
modulo a re-scaling of the volatility of volatility. The first-order asymptotic is used in Corollary 13.161 below to 
highlight differences with the Heston model. 

Let us now define the following functions: 
(3.19) 

Mt, r) := T (^^f (p2 + 2) + tidV^ + ^v{^p - in)^ + 2Kt^^ (O - V^) , 
Pir (e^l - - 8kw + 2^pv) fpt 



1^2 (i, r) ( (521p^ - 452p2 + 56) f + mK0y/^ {2p^ - l) + 40p\(/3C - 4k) 

f ('i /I „2 r-^ ,-2 , n.^n /-/in ■7^2\ , /'7„2 a\ \ , ^^^^ 



7680.2 



48„2 ( (1 V -5)e + 2neV^ (10 - 7p') + 2nv {7p' - 4) ) + A_ + 4^7^ {V^ ~ 26 

In order to gain some intuition on the role of the Schobel-Zhu parameters on the forward smile we expand 
our solution (to first order in e) around the ATM point in terms of the log strike k. 

Corollary 3.16. The following expansion holds for the Schobel-Zhu forward smile as e and k tend to zero: 
(3.20) 4,,,(fc) = v + eMt,T) + (^-l + eu,{t,T)^ k+ (^^^^^ + ^+e,,,{t,T)^ ^ 

+ 0ik^) + 0iek^) + 0{e^), 
where vq, vi and V2 cii^e defined in p.lQp . 

Remark 3.17. At this order we can make the following remarks concerning the SZ forward smile: 

(i) Remark l3.4r iv) for the Heston forward smile also applies to the SZ forward smile. 

(ii) The forward ATM volatility has a different dependence on the volatility of volatility ^ in Heston and SZ. 
In Heston (Remark l3.4r u)). a-^t.eriO) — cro,eT(0) is decreasing in ^. In the SZ model, Corollarv 13.161 implies 

<7et,eT{0) — o'o,er(0) = (9 ~ \/v)Kte + 0{e'^), as e tends to zero. 
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which does not depend on ^ (up to an error of order 0{e^)). Also for realistic parameter choices (eg. 
p < 0) the Heston ATM forward volatility is decreasing in ^ while in the SZ model (for example when 
£, > 2v) it is increasing in ^ and the impact is small. This effect is illustrated in Figure [31 

Remark 3.18. An analysis analogous to that of the Heston model can be conducted for the large-maturity SZ 
forward smile. We shall omit it here though for brevity. 





(a) Heston forward smile ^ dependence 



(b) SZ forward smile ^ dependence 



Figure 3. Here t = 1/2 and r = 1/12 and wc apply Corollaries 13.161 and 13.31 Circles use the 
Heston parameters v = 0.07, 9 — 0.07, k — 1, p — —0.6, ^ = 0.3 and SZ parameters v = 0.07, 9 = 
\/0.07, K = I, p = —0.6, ^ = 0.3. Squares use the same parameters but with ^ = 0.4. 



3.4. Time-changed exponential Levy. It is well-known (see for example in [T4j Proposition 11.2]) that the 
forward smile in exponential Levy models is time-homogeneous in the sense that at^r does not depend on t for any 
fixed r > 0, since the process has stationary increments. This is not necessarily true in time-changed exponential 
Levy models as we shall now see. Let (Yt)t>o be a Levy process with mgf given by logE (e"^') = t(l){u) for 
all i > and u £ T)^ := {u eM. : |0(u)| < oo}. Wc consider models where (Xt)t>o {Yvt)t>o pathwise and 
the time-change is given by Vt := /q Vsds with v being a strictly positive process independent of Y . We shall 
consider the two following examples: 



(3.21) 
(3.22) 



dvt ~ —Xvtdt + dZi, 



with vq ^ V > and k,^,9,X > 0. Here VF is a standard Brownian motion and Z is a compound Poisson 
subordinator with exponential jump size distribution and Levy exponent l{u) := Xdu/{a — u) for all w < a 
with d > and a > 0. In (|3.2ip . v is a Feller diffusion and in p.22|) . it is a F-OU process. We now define the 
functions V and H from Poo to R by 

.,^^(,^^;;t^2^\ and 



(3.23) V{u) 



and H{u) 



2^9 



k9 - 2PtV{u) e 



los 



k9 - 2PtV{u) 



K0(l-7(0(w))) 



and the functions V and H from D^o to 

0(u)Ad 

V (U) := - 

(3.24) 



by 



Viu) 
Hiu) 



a\ — 4'{u) ' 
Xad 



log 1 - 



4>{u)\ (l){u)ve 



-At 



A 



dlog 



i:Ae 



At 



atA 



((/)(u) — aA) 
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where we set 

(3.25) := {u : dp{u) < AiV(2e')} , := {u : dp{u) < a\} , 

(j) is the Levy exponent of Y and j3t and 7 are defined in (j5.46p . The following proposition — proved in Sec- 
tion 15.51 — is the main result of the section. 



Proposition 3.19. Suppose that (p is essentially smooth fDefinition \2. 2\) . strictly convex and infinitely 
entiable on with {0, 1} C I?^ and 0(1) = 0. Then Corollary \2.10\ and Proposition \2. l4\ hold: 

(i) if V follows p.21[) . with Aq.o ~ V , Ai.o = H, A2,o ~ and Vq^q ~ Poo/ 

(ii) if V follows p.22p . with Aq^o = V, Ai^o = H, A2,o ~ and I?o,o = T^oc! 
(Hi) if vt = 1, with Ao,o ~ 0? ^1.0 = 0, A2.0 = and I'o.o = '^(p- 

Remark 3.20. 

(i) The uncorrelated Heston model p.ip can be represented as Yt := —t/2 + Wt time-changed by an integrated 
Feller diffusion (|3.2ip . With (j){u) = u{u—l)/2 and I?^ = M, Proposition l3.19r i) agrees with Proposition l3.8l 

(ii) The zeroth order large-maturity forward smile is the same as its corresponding zeroth order large-maturity 
spot smile and differences only emerge at first order. It seems plausible that this will always hold if there 
exists a stationary distribution for v and if v is independent of the Levy process Y; 

(iii) Case (iii) in the proposition corresponds to the standard exponential Levy case (without time-change). 

We now use Proposition l3.19l to highlight the first-order differences in the large-maturity forward smile (|2.17p 
and the corresponding spot smile. If v follows p.2ip then a Taylor expansion of vf' in (|2.16p around t = {) gives 

vT{t,k)-vT{^. k) = J^llf^j^y Vium [^^^^^^ + 1 - ^) t + 0{t% for aU k e M\ {7' (0), 7' (1)}. 

Using simple properties of ug° discussed below ()2.16p and V we see that the large-maturity forward smile is 
lower than the corresponding spot smile for k G {V (0), V (1)) (which always include the at-the-money) iiv>6. 
The forward smile is higher than the corresponding spot smile for k S M\(y'(0), ^^'(1)) (OTM options) if u < 0, 
and these differences are increasing in ^/k and t. This effect is illustrated in Figure |4] and k G {v' (0) , v' {1)) 
corresponds to strikes in the region (0.98, 1.02) in the figure. 

If V follows ()3.22p then a simple Taylor expansion of v^{-, k) in (|2.16p around t = gives 

Similarly we deduce that the large-maturity forward smile is lower than the corresponding spot smile for k G 
(y (0), V'{1)) if u > d/a. The forward smile is higher than the corresponding spot smile for k € M\(y (0), V'{1)) 
(OTM options) if w < d/a, and these differences are increasing in t. 

If V follows p.2ip f p.22p ) then the stationary distribution is a gamma distribution with mean 9 {d/a), 
see [m page 475 and page 487]. The above results seem to indicate that the differences in level between the 
large-maturity forward smile and the corresponding spot smile depend on the relative values of vq and the 
mean of the stationary distribution of the process v. This is also similar to Remark I3.4f ii) and the analysis 
below Remark 13.91 for the Heston forward smile. These observations are also independent of the choice of </> 
indicating that the fundamental quantity driving the non-stationarity of the large-maturity forward smile over 
the corresponding spot implied volatility smile is the choice of time-change. 
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As an example of a Levy process satisfying the assumptions of Proposition l3.19l consider the Variance- Gamma 
model, defined by 



(it) = /iu + Clog 



GM 



for all u e (-G,M), 



{M ^ u) {G + u) J ' 

with C > 0, G > 0, M > 1 and fi — Glog ^ (m-Txc+i) ) "insures that (e^')t>o is a true martingale (0(1) = 0). 
We immediately obtain 

G(G - M + 2m) + ^(G + m)(M - u) 



and 



0"(«) = 



G ((G + M)2 + {2u + G- Mf) 



{G + u){M~u) 2(G + u)2(A/-u)2 

for u G (— G, M), so that is essentially smooth and strictly convex on (— G, M). It is also infinitely differentiable 
on (— G, A/) with {0,1} C {~G,M) and Proposition 13.191 applies. For Proposition I3.19r iii) we can compute 
u* : M ^ (-G, M) through in closed-form. The solutions to (t>'{u*{k)) = k are m*(/^) = (A/ - G)/2 and 



< (fc) ^ - - - ^) % V^^f ^ + + - for all 



The sign condition (A/ - u) {G + u) > imposes -2G± ^4G2 + {G + AI)^ {k - n)^ > for all fc 7^ pt. Hence 
(continuous on the whole real line) is the only valid solution and the rate function is then given in closed-form 
as A*(fc) ~ ku\[k) — (f){u\{k)) for aU real k. 





0.8 0.9 1.0 1.1 1.2 1.3 0.8 0.9 1.0 I.I 1.2 1.3 

(a) Feller time-change: forward smile vs spot smile u > 8. (b) Feller time-change: forward smile vs spot smile d < d. 

Figure 4. Circles represent t ~ Q and t = 2 and squares represent t = 1/2 and t = 2 using 
a Variance-Gamma model time-changed by a Feller diffusion and the asymptotic in Proposi- 
tion [nH In (a) the parameters are G = 58.12, G = 50.5, M = 69.37, k = 1.23, 6 ^ 0.9, 
^ = 1.6, u = 1 and (b) uses the same parameters but with 6 = 1.1. 



4. Numerics 

We compare here the true forward smile in various models and the asymptotics developed in Propositions l2.1^ 
and 12.141 We calculate forward-start option prices using the inverse Fourier transform representation in [41[ 
Theorem 5.1] and a global adaptive Gauss-Kronrod quadrature scheme. We then compute the forward smile at.r 
and compare it to the zeroth, first and second order asymptotics given in Propositions 12.1^ and \^I7lM for various 
models. In Figure [5] we compare the Hcston diagonal small-maturity asymptotic in Proposition 13.21 with the 
true forward smile. Figure [5] tests the accuracy of the Heston large-maturity asymptotic from Proposition 13.81 
In order to use this proposition we require < p < min (p_|-, k/^) with p± defined in p. lip . For the parameter 
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choice in the figure we have p_ = —0.65 and the condition is satisfied. Finally in Figure [7] we consider the 
Variance Gamma model time-changed by a F-OU process using Proposition 13.191 Results are in line with 
expectations and the higher the order of the asymptotic the closer we match the true forward smile. 



FwdSmile 




Strike 



(a) Heston diagonal small-maturity vs Fourier inversion. (b) Errors 

Figure 5. In (a) circles, squares and diamonds represent the zeroth, first and second order 
asymptotics respectively in Proposition 13.21 and triangles represent the true forward smile us- 
ing Fourier inversion. In (b) we plot the differences between the true forward smile and the 
asymptotic. We use t = 1/2 and t = 1/12 and the Heston parameters v = 0.07, 9 = 0.07, 
K= 1, ^ = 0.34, p= -0.8. 



FwdSmile 

Error 




Striiie 



(a) Heston Large-Maturity vs Fourier Inversion. (b) Errors 

Figure 6. In (a) circles, squares and diamonds represent the zeroth, first and second order 
asymptotics respectively in Proposition 13.81 and triangles represent the true forward smile us- 
ing Fourier inversion. In (b) we plot the differences between the true forward smile and the 
asymptotic. We use t — 1 and t = 5 and the Heston parameters v — 0.07, 6 — 0.07, k = 1.5, 
^ = 0.34, p = -0.25. 
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(a) F-OU time-change large-maturity vs Fourier inver- 
sion. 



Error 
0.005 r 





(b) Errors 



Figure 7. In (a) circles, squares and diamonds represent the zeroth, first and second order 
asymptotics respectively in proposition 13 . 1 9l and triangles represent the true forward smile using 
Fourier inversion for a variance gamma model time-changed by a F-OU process. In (b) we plot 
the differences between the true forward smile and the asymptotic. We use t ~ 1 and t — 3 
with the parameters C = 6.5, G = 11.1, AI = 33.4, v ^ 1, a = 0.6, d = 0.6, A = 1.8. 



5. Proofs 

5.1. Proofs of Section [2l We define a change of measure by 



(5.1) 



dP 



exp 



with u*{k) defined in By Lemma [llUi), u*{k) e 2?go for all fc G R. Since Pg^ C 2?^ this means that 

Ae (ii*(fc)) < oo. Also dQfc_e/dP is almost surely strictly positive and by the very definition E(dQA;/dP) = 1. 
Therefore ()5.ip is a valid measure change for all fc e R. Wc define the random variable 

(5.2) Zfc,, {X, - k) l^e 

and set the characteristic function ^ : M C of Z^.e in the Q^.e-measure as follows 

(5.3) $z,,,(u) =E^'='^ (e'"^'=--) . 

Recall from Section [2] that for ease of exposition A,.; := A^^j (w*(fc)) with A.^ ; defined in (|2.4p . This notation 
will be used throughout the section. Wc now have the following important technical lemma. 



Lemma 5.1. The following expansion holds for the characteristic function e defined in 
log$z..,(u) = -^Ao,2u^ + m{u)Ve + V2{u)e + iMe^^^ + O (e^) , 



as £ — > 0; 



where the functions rji, i = 1, 2, 3 are defined in (|5.4p . 



Remark 5.2. By Levy's Convergence Theorem [SH Page 185, Theorem 18.1], Zk.e defined in (|5.2p converges 
weakly to a normal random variable with mean and variance Ao,2 in the Q^^g-measure as e tends to zero. 
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Proof. By Lemma [2^ u*(fc) G 'Dq q, therefore using the definition of the Qfe ^-measure in (|5.ip we have 



log$z,,(u) =logEF 



logE^ 



cxp 



cxp luye 



iA,K(fc))-^+logE' 



cxp 



luve + u 



(fc)) 



-^ + i(A,(iuVi + M*(fc)) - A, {u*{k))) 

ik 1 



= — ^ + - XI (^0'" + '^^i'" + '^^'^2,n + O (e^)) (iux/e)" , 

^ 71 — 1 

where the last equahty holds by the differentiation and expansion properties in Assumption 12.11 We now write 

\n—Z n—1 



log$z,,,(M) 



ilifc iAn ^U 1 . 9 1 

^ + — ^--AoV + - 



= 4w + E^( 

n— 1 ^ 



,(n + 2)! 



where we used (12.51) from Lemma 12.41 and where 



(5.4) 



7?i(u) iwAi^i - iiw^Ao,3, 

1 1 

miu) -^w^Ai,2 + — u'^Ao,4, 

7?3(m) := iuA2,i - -iM^Ai_3 + — iu'^Ao,5- 



Proof of Theorem \2.6\ For j = 1, 2, 3, let us define the functions gj : by 

Si(2;>y):="( if.? =2, 

min(a;, y), if j = 3. 

Using the definition of the Q^^j-measure in (jS.ip we can re-write the option price as 



-9j 



C7Ae(«'(fc))]gQfc, 

= Q-jiku' {k)-A,{u- (k))]^Q^,, 

By the rescaled mgf expansion in Assumption 12.11 and Equality (j2.6p we immediately have 
(5.5) exp (^-i(fcu*(/c) - Ae(u*(fc)))^ = cxp (^-^A* (fc) + Ai,o + Az.o £ + O (e^)^ 



From the definition of the random variable Zh ^ in (15.21) wc have 



where for j = 1, 2, 3, we define the modified payoff functions : R — > M+ by 



□ 
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For a function / G i^(M) denote its Fourier transform by [F f){u) := /jj e^"^/(z)dz, for any w e R. Assuming 
(for now) that gj E L^(R), we have for any u G R, 



" exp l'-^^^^) g, (e'^f^'\l) e^"Mz, for j = 1, 2, 3. 



For j = 1 we can write 

exp (^z {^/ef (e) ~ u*{k)/y/s + ±u) 

{u*{k) ~ ±u^/e) (ii*(fc) — ef{e) — ±u^) ' 

which is vaUd for u*{k) > ef (e). For e sufficiently small and by the definition of / in (|2.7p this holds for 
u*{k) > c. For J = 2 we can write 



exp (^z {—u* (fc) / ^/e + iu) ^ 
—u*{k)/y/e + ±u 



exp (^z{-u*{k)/^/e + iu) 
~u*{k)/^/e + iu 



exp 



{e)-u*{k)/^+iu)) 
y/ef (e) - u*{k)/y/e + iu 



{u*{k) — iuy/e) (u*(fc) — £./(£) — iuy/e)' 
which is valid for u* (fc) < as e tends to zero. Finally, for = 3 we have 







g3(2)e'"Mz= / e-'^'^gs (e'^^(^\l) e^"^dz + J e'^'g^ (e'-^f^'^l) e^-^'d 





oo 



Z 



exp (^z {y/ef (e) - u* {k)/^/e + iu) 
v^/ (e) - It* {k)/y/e + iu 

e'/'m 

(ii*(fc) - iuy/e) {u*{k) - ef{e) - iu^/e) ' 



exp (^z {-u* (k) / y/e + iu)^ 

—u*{k)/y/e + iu 



which is valid for < u*{k) < ef (e). For e sufficiently small and by the assumption on / in ()2.7p this is true for 
< u*{k) < c. In this context u*{k) comes out naturally in the analysis as a classical dampening factor. Note 
that in order for these strips of regularity to exist we require that {0, c} C Vq q, as assumed in the theorem. By 
the differentiability property in Assumption 12.11 and the strict convexity and essential smoothness property in 
Assumption 12.31 we have 



(5.6) 



0<u*{k)<c if and only if Ao.i(O) < fc < Ao,i(c), 
u*{k)<0 if and only if fc < Ao.i(O), 
u*{k)>c if and only if fc > Ao.i(c). 



For e sufficiently small and the strips of regularity defined above, the modified payoffs gj are in i^(R). By 
Remark l5.2[ Z^.e converges weakly to a zero-mean Gaussian random variable as e tends to zero and the Gaussian 
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density is in i^(R). For e sufficiently small we can apply Parseval's Theorem [29l Theorem 13E] to write 



(5.7) EQ*-- MZk,e)] = <^ 



2tt (u*(fc) + iu^/e) {u*{k) - ef{e) + iu^/e) 
1 f e'/^f{e)<fz,Au)du 



2-K 7k (ii*(fc) + luy/e) {u*{k) - ef{e) + luy/e) 
1 /■ e^'^f{e)^z,Au)^u 



. 27r 7k (u*(A:) + ±u^/e) {u*{k) - e/(£) + ±Uy/e) 



if j = 1, ?/*(/(:) > c, 
, if J =2,w*(fc) <0, 

if i = 3,0 < u*{k) < c, 



where we have used that 



e3/2/(£) 



{u*{k) ~ ±uy/e) — e/(e) — iuA/e) (u*(fc) + ±uy/e) {u*{k) — ef{e) + iuy^) ' 

with a denoting the complex conjugate for a G C. Now using Lemma |5. II we write 

£3/2/(e)$(if)dM f f AoW" 



(5.8) 



/ 



{u*{k) + ±uy/e) (u*(fc) — ef{e) + ±uy/e) 
where the function H : M_|_ x M — >■ C is defined as 



exp 



H (e, u) du, 



_ exp {m{u)V^ + T]2{u)s + m{u)e^^^ + O (e^)) x 3/2 
■ (u* (fc) + iw^i) (u* (fc) - efie) + ±uV^) ^ ^^'^ ' 

with rji (i = 1,2, 3) defined in (|5.4p . A Taylor expansion of i? around e = for c = yields 



H 



3/2 



U*(fc)2 

where we define the following functions: 
(5.9) 

hi{u,c) 

hi{u,c) : 



I + hi{u,0)y^ +h2{u,0)e + h3{u,0)s^^^ + + ( ^ 



(fc) \u*{k) u*(kf 



*{k) - c \u*{k) 



/i2(u,c) 



(c2 - Zcu*{k)+iu*{kY) 
u*{kY{u*{k)-cf 



iw3(4it*(fc)3_c3 + 4c2w*(fc)-6cu*(fc)2) . , , , , ±u [c - 2u* {k)) 



h2{u,c) 
hs.{u,c) 



u*{kf {u*{k)~cy 

-^^^+??2(w) ^ 9 ^ ' 

2 u*{kY{u*{k)^cf 
— h 7Ji(U)r]2(u) + 773 (u) + 



u*{k) — c \w*(fc) 
n ,,3 



u*{k) (u*(fc) 



- 2 



itt / V^M , iuriiUi) 



*{k)-c 



(fc)3 u*{kY {u*{k)-cf 

Iv!^ ~ v^u* {k)rji{u) ±u f r]i[u) 



(fc) u*{kYJ u*{k){u*{k)-cf u*{k)\ 2 



'72 (w) 
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with the rji for i = 1, 2, 3, defined in (|5.4p . A Taylor expansion of H around e = for c > yields 
Hie,u) = 

{l + Vi{u)V-e + + V2{u)^ e + + V2{u)vi{n) + e'/' + O (e^) } 

'(fc)(£/(£)-e) 2±u^{ef{e)-c)'. 
c{u*{k)-c) c{u*{k)-cf 



u*{k) {u*{k) - c 



, / ^ ^ , , X , / \ ■1/9 u*(k)(ef(e) - c) 2±uJe{e f ie) - c)u* ik) , 



u* [k) [u* [k] — c) 1^ ~'~ 1^ ' / ^ ' ' ^ ' ■'^ ' / ' c{u*{k) — c] 

with the ft-i, ft-i defined in (|5.9p and the r]i defined in (j5.4p . We will shortly be integrating H against a zero- 
mean Gaussian characteristic function over ffi. and as such all odd powers of u will have a null contribution. In 
particular we note that the polynomials 

and 

\u*ik) {u*{k)f j ^' c{u*ik)-c) V u*{k)~cj 

are odd functions of u and hence have zero contribution. The major quantity is /i2, which we can rewrite as 

(5.10) h2{u,c) = h2,i{c)u + h2.2{c)u 

where 

2Ai4 cAi,i Af i+Ai.a - 3cit*(fc) + 3u*(/fc)2 



/i2,i(c) := 
/i2,2(c) := 



Let 



W*(fc)-C U*(fc)(w*(fc)-C) 2 M*(fc)2(M*(fc)-c)^ ' 

cAo^3 Ao^3 Ai^iAp^a Ao^4 

6M*(fc)(M*(fc)-c) 3{u*{k)-c) 6 24' 

cVil{c>o} + e3/V(e)l{c=o} 



- c) 

Using simple properties of moments of a Gaussian random variable we finally compute the following 



/ 



exp I H (e, u) au 



In the second line we have dropped all odd powers of u and have used (|5.10p . Combining this with (|5.8p . ()5.7p 
and (|5.5p with the property ()5.6p . the proposition follows. □ 



In |25| Gao and Lee have obtained representations for asymptotic implied volatility for small and large- 
maturity regimes in terms of the assumed asymptotic behaviour of certain unspecified option prices, outlining 
the general procedure for transforming option price asymptotics into implied volatility asymptotics. 

The same methodology can be followed to transform our forward-start option asymptotics (Corollary 12.81 
and Corollary I2.10p into forward smile asymptotics. In the proofs of Proposition 12.121 and Proposition 12.141 
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we therefore assume for brevity the existence of an ansatz for the forward smile asymptotic and solve for the 
coefficients. We refer the reader to [IS] for the complete methodology. 

Proof of Provosition \2.12l Using Ao,i(0) ~ and substituting the ansatz ^^(^^(fc) = vo{k,t,T) +ui(fc,i,r)e + 
V2{k, t, r)£^ + O (e^) into Corollarv 12.91 we get that forward-start option prices have the asymptotics 



E 



( e " 



Hk>0} 



(e^-e^i^^")" 



Hk<o} 



cxp 



fc2 



Pvi{k,t,T) k\ {vo{k,t,T)eT)^^^ 



k^V2^ 



{l+j{k,t,T)e + 0{e')) 



where we set 



l{k,t,T) := ~TVoik,t,T) ( 77 + ^ 1 + 



k'^V2{k,t,T) k'^vi{k,t,TY 2>vi{k,t,T) 



+ 



>2 sy 2TWo(fc,i,T)2 2r?;o(fc,i,r)3 2wo(fc,<,T)' 
The result follows after equating orders with the general formula in Corollarv 12.81 



□ 



Proof of Proposition \2.14\ Substituting the ansatz 

(5.11) al,{k)=v^{k,t)+vrik,t)/T + v^{k,t)/T^+0{l/T^), 

into Corollarv 12 . 1 1 1 we obtain the following asymptotic expansions for forward-start options: 



E 



(e^^"-e^- 
-r( 



exp 



\2vo{k,t) 
4T-'/h,o{k,tf/^ 
(4fc2-l;o(fc,^)2)^/2^ 



1a -E 

k 
2 



e^^" A e^^ 



1b + E 



Ic 



vo{k,t)\ vi{k,t)k^ vi{k,t) 



8 



2vQ{k,tY 



1 



8 



where 
(5.12) 



A 



7-(fc,t) 



C:=|fc<-ia,%(fc)|, 
I2k^ + vl{k,t)){Ak'^vi{k,t)^vl{k,t){vi{k,t)+^)) v^ik,t)P V2ik,t)k^ V2{k,t) 



2vo{k,t){vl{k,t)-4k^f Hik,t) 2v^{k,t) 8 

We obtain the expressions for uj" and by equating orders with the formula in Corollarv l2.10l However it is 
not clear which is the correct root for the zcroth order term . In order to do so, we have to match the domains 
in (j5.12p and in Corollarv 12.101 Indeed, suppose that we choose the roots according to in ()2.16p . For t 
sufficiently large the condition k > a'^/2 is equivalent to fc > wg°(fc,<)/2. Now for k > Ao,i(l) or fc < Ao,i(0), 
the definition of in (|2.16p implies 



(5.13) 



fc > cr 2/2 if and only if ^J {A* (k) - kf + k {A* (k) - k) > A* (fc) - fc. 



which is always true since A*(fc) > fc by Lemma [2.4f iii). Now, for fc G (Ao_i(0), Ao,i(l)), the definition of 
in (j2.16p implies 



(5.14) 



fc >crV2 if and only if - [A* {k) - kf + k {A* [k) - k) > A* (fc) - fc. 



which never holds. By the assumption in the proposition 12.141 and Assumption 12.11 we have {0, 1} C T>q q and 
Ao.o(O) = Ao,o(l) — 0. The differentiability and strict convexity of Aq.o (Assumptions [2731 and [2A|) then imply 
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Ao,i(0) < and Ao,i(l) > 0. Since > we can ignore the case k < Ao,i(0) < and hence k > if and 
only if fc > Ao.i(l). Similarly the definition of in (|2.16p implies that for t large enough r, 

< /c < if and only if Ao4(0) < fc < Ao,i(l) and fc < -ct^/2 if and only if fc < Ao,i(0). 

This lines up the domains in ()5.12p with the domains in Corollarv l2.10l Had we specified the roots in any other 
way, it is easy to check that a contradiction would have occurred. □ 

5.2. Proofs of Section 13. IL For a stochastic process {Xt)t>o we define the forward mgf of the process by 

(5.15) A(u) := logE fe"^-*') , for all u e Pa, 



where 'D\ := {u € R : |A(u)| < oo} and Xr '^ is defined in (|2.1ip . We now let {Xt)t>o be the Heston process 
satisfying the SDE p.ip . By a straightforward application of the tower property for expectations (see also [51]). 
the forward mgf defined in (|5.15p is given by 

B (u t) 2k9 

(5.16) A{u)^A{u,t) + - ofl p/ . fe-"*- — log(l-2/3t^(^,r)), for all u € I^a, 

where 

1 — "f{u) cxp (— (i(u)T) 

^ I": ' ; •= — I l'^ ^ Z^?'" - - ^i<Jg I 

(5.17) 



B{u,t) 



K — p^u — d{u) 1 — exp (~d(u)T) 



1 — 7(m) exp (— (i(u)r) ' 

and d, 7 and /3 were introduced in p.l3p . In the next two sections we develop the tools needed to apply 
Propositions 12.121 and [2TT41 to the Heston model. 

5.2.1. Proofs of Section \3.1.1\ Wc consider here the Heston diagonal small-maturity process [X'frA with X 

V / £>0 

defined in ()3.ip and [xi^^] in (|2.1ip . The forward rescaled mgf A^ in \2.1\ is easily determined from (|5.16p . 



r>0 

Lemma 5.3. For the Heston diagonal small-maturity process we have I?o,o ~ T^t.r o,nd G 2?qo with "Dt^r 
defined in (|3.3p and I?o,o fw^'^ ^00 defined in Assumvtion \2.1\ 



Proof. For any i > 0, the random variable Vt in (|3.ip is distributed as j3t times a non-central chi-squarc random 
variable with q = 4k6'/^^ > degrees of freedom and non-centrality parameter A = vc~'^*/(3t > 0. It follows 
that the corresponding mgf is given by 

1 



, for all u < 



(5.18) Ar(.):^E(e-)^exp(^)(l-2A.)- , 

The re-normalised Heston forward mgf is then computed as 



E 



QTi^et + ET—^Et) 



(e^(^-+"-^-)|J-,,)] = E (e^(f ^-)+^(f ^-)^-) = e^(t-^^)A,^, {B {u/e^er)) , 



e 

which agrees with (|5.16p . This only makes sense in some effective domain 'D^t,eT C M. The mgf for Vet 
is well-defined in := {u G M : B{u/e,eT) < and hence Vet,er = n V^^, where V^^ is the 

effective domain of the (spot) Heston mgf. Let us first consider 2?^ for small e. From [Sj Proposition 3.1] if 
^^(u/e — l)u > {k — S^pu/eY then the explosion time t*{u) := sup{t > : E(e"'^*) < 00} of the Heston mgf is 



(7) ^ Jf^(u/.^^^ul^(.^nfu/.Y2 -l{PCVe-K<o} +arctan 



^e{u/e-l)u/e-{K-p^u/ey 



^e{u/e~l)u/e^{K~ piu/eY \ ^"^^'^ '^^'■'^ I p^u/e-K 
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Recall the following Taylor series expansions, for x close to zero: 



arctan ( — - — y \l f — — l) — — (k — CP") ) = sgn(u) arctan [ - ) + O (x) , if p 7^ 0, 



P^u/x — k\ \x / X \ X 



aictan (-^J e ~ ^ - nA =-^ + 0{x), if p = 0. 




As £ tends to zero £,'^{u/e — l)u/e > (k — p^u/e)^ is satisfied since > ^^p^ and hence 

rh {u/e) = { £.\u\ V''^{''=°> + I {^hpum + sgn(u) arctan ^ l{p_,o} + 0{e)^ , if u ^ 0, 

if M = 0. 

Therefore, for e small enough, we have (^) > er for all m € (u-,^^.), where 

u_ ^ arctan l{p<o} - ^l{p=o} + ^ (^arctan 0^ - tt^ l{p>o}, 

w+ := ^ (^arctan (^Q + tt^ l{p<o} + ^l{p=o} + ^ arctan l{p>o}. 

So as £ tends to zero, 2?|^ shrinks to Regarding 'D^^, wc have (see (|5.22p for details on the expansion 

computation) P^tB [u/ e^er) ~ ^A(it,0, r) + 0(e) for any u € (u_,u+), with A defined in (|3.2p . Therefore 
lim^N^o^^rt = e R : A(m, 0,t) < ^} and hence \in\\^o V^t.er = {u £ R : A (w, 0, r) < ^} n (m_,m+). It 
is easily checked that A {u, 0, r) is strictly positive except at u ~ where it is zero. A' (u, 0, r) > for u > 0, 
A' (u,0, r) < for u < and that A(u,0, r) tends to infinity as u approaches u±. Since v and ^ are strictly 
positive and i > it follows that {w G M : A(u,0, r) < ^} C (u_,w+) with equality only if t = 0. So 2?o,o is 
an open interval around zero and the lemma follows with I?o,o — T^t.r- D 

Remark 5.4. For w G M* the inequality < A(w, 0,t) < ^ is equivalent to A(it,t,T) S (0,oo). where A is 
defined in (|3.2p . In Lemma TS.SI below we show that A is the limiting mgf of the rescaled Heston forward mgf 
and so the condition for the limiting forward domain is equivalent to ensuring that the limiting forward mgf 
does not blow up and is strictly positive except at u = where it is zero. 

Lemma 5.5. For any t > 0, t > 0, u G Vt r, the following expansion holds as e tends to zero: 

K,{u) = A(u, i, t) + L{u, t, T)e + O (e^) , 
where Vt.T, A and L are defined in (|3.3p . p.2p and ^AA\ and A^ is the rescaled mgf in Assumvtion \2. 1\ for the 



Heston diagonal small-maturity process \X^'^r\ 

V /e>0 



■ e>C 

Remark 5.6. For any u G I?t,T! Lemma 15.31 implies that h.^{u) is a real number for any e > 0. Therefore L 
defined in p.4p and used in Lemma 15.51 is a real- valued function onT^t t- 



Proof. All expansions below for d, 7 and fit defined in p.l3p hold for any u G Vt^r'- 

d {u/e) = - (n^e^ + ue - 2kp) - u^^^p) = —do + di+0{e), 
e e 



= .e-p^u + ±udo + d,e + 0{e-) = --9.+0{e) 



(5.19) /3^, ^ i^2^£ _ + O (e^) 
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where we have set 

^r; on^ ^ -c ^ ^ ^ i {2np - Q sgn(M) ip- psgn(M) (2k - gp) sgn(M) 

5.20 do pCsgn w , rfi — , go — gi — 

2p ip + psgn(w) ^p(p+ ipsgn(ii))^ 

where sgn(w) = 1 if m > and -1 otherwise. From the definition of A in (j5.17p we obtain 
(5.21) A (w/£, et) = — f (k - pe^/e - d(u/e)) er - 2 log f V- 7Kg) 



— ((ke - p^M - itido - edi + ©(e^)) t 



-2 log 



1 ~ (ffo — ^£gi/u + 0(£^)) exp (— iudoT — edir + 0(e^)) 



l-{go-iegi/u + 0{e^)) J J 

= Loiu,T)+0{e), 

where Lq is defined in (j3.5p . Substituting the asymptotics for d and 7 above we further obtain 
1 - exp (-d(u/e)er) 1 - exp (-iudoT - £diT + 0{e^)) 



1 — j{u/e) exp {—d{u/e)eT) 1 — (50 — + C(e^)) exp (— iudgr — edir + 0{e'^)) ' 

and therefore using the definition of B in (|5.17p we obtain 

K ~ p^u/e — d{u/e) 1 — exp (— d (w/e) er) 



(5.22) B{u/e,eT) 



1 — 7 (u/e) exp (— d (w/e) er) 
p£^u + ±udo 1 — exp (— iudoT) 



1 ^ .90 exp (— iudor) 

A(u, 0,r) ^ , , . 

+ Li(u,r) +©(£), 



+ Li(u,T)+0(e) 



we 

with Li defined in ((3?5|) and A in ((3?2l) . Combining (|5.19p and (|5.22p we deduce 

(5.23) j3,tB{u/e,£T) ^ — +( — le + 0(e 

and therefore as e tends to zero, 

eB(M/£,£T)we-'=^* _ [A(u, 0, r) + wLi(w, r)e + O (e^)] (l - t< + ©(e^)) 



1 - 2l3etB{u/e, er) 1 - S,HK{il, 0, t)/2w + (A(u, 0, T)KS,H'^/'iv - Li{u, T)^^t/2) e + O (e^) 

/ vLi{u,t) K^'^t- 
VA(m,0,t)2 4^ 



(5.24) = A(7/, t,T)+[ Aiu, t, rf f ^^^^^ - ^ ) - AitA(?/, t, r) ) e + ©(e^). 



Again using ()5.23p we have 

(5.25) log (1 - 2/3etS (u/e, er)) = -—\og[\~ ^ j £ + ^(e'). 

Recalling that 

A, u = eA w/e, er + ' , ' . ve - —- log 1 - 2/3^*5 w/e, er , 

the lemma follows by combining (|5.2ip . (|5.24p and (|5.25p . □ 
Lemma 5.7. For aUt>0,T>0,A (given in (|3.2p j ?s convex and essentially smooth on Vt^r, defined in (|3.3p . 
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Proof. The first derivative of A is given by 



p( + ^C^tu- cot {^^pTu) {-(_p - l^Hu + ^p cot {^^pTu)y 

K(u,t,T) ^ A{u,t,T) ( ^ 1 2-2_ 2/^1 



du 

^ ^" ^ \ + ( i_ + i^^pV csc^ ( -p^Tu 

Any sequence tending to the boundary satisfies A(u,0, t) — > 2vl^t which imphes A(M,t,T) /• 00 from Re- 
mark [5]4] and hence |9A(w,t,T)/9w| 00. Therefore A(-,t,r) is essential smooth. Now, 

a2A(w,t,T) _ 2ww(~^-ie'Acsc2(V'0) e'i^, + p2^Vv(l-^„cot(7/>„))csc2(V'„) 



{-ip - \etu + pecot(V'„))' ■ {-pi - ^etu + pecot(V'„))' 

+ pVcsc2(V'„))^ i; + pVi;(l-'0HCot(V'„))csc2(V'„) 

= —A(U,t,T) 2 2 ■ 

(p + - pcot{ipu)) {p + ^(.tu - pcot(V'„)) 

where ipu ■= \p£,Tu. For u e 2?i,T- \ {0}, we have K{u,t,T) > and A{0, t,T) = from Remark 15.41 Also we 
have the identity that 1 — 6/2 cot {0/2) > for 6* € {—2tt, 2tt), so that A is strictly convex on T>t,T- O 

Proof of Corollarv \3.3[ Wc first look for a Taylor expansion of u*{k) around k = using A' {u* {k),t,T) = k. 
Differentiating this equation iteratively and setting = (and using u* (0) = 0) gives an expansion for u* in 
terms of the derivatives of A. In particular, A"(0, r)u*' (0) = 1 and A"'(0,t,r)(u*'(0))2 + A"(0,i,r)u*"(0) = 0, 
which implies that u*'(0) = l/A"(0,t,T) and u*"(0) = -A"'(0, i, t)/A"(0, r)^. From the explicit expression 
of A in (|3.2p . wc then obtain 

u*(k) = A _ ^fc2 , CM(lV-4)r-12t) ^^3 5e^p(48t+(l6-37p^)r) ^^, 
^ ' TV 4tv^ 24t2w3 192t2w4 

(1080^2 + (2437p^ - 1604p2 + 112) ~ 180 (27p^ ~ 4) rt) 

1920r3i;5 ^ ''■ 

Using this series expansion and the fact that A*(fc) = u* {k)k — A(it*(fc), t, r), the corollary follows from tedious 
but straightforward Taylor expansions of Vq and Vi defined in (|2.14p . □ 

Lemma 5.8. Under the stopped- share-price measure (j2.18p the forward Heston mgf defined in (|5.15p reads 

A{u) = A{u, t) + ^1^2^! „c-«* - ^ log f 1 - 2%B{u, t) \ , for all u G I?a, 

i-2ptB{u,T) r ^ ^ 

where A and B are defined in (|5.17p . 

~ ^2 

(5.26) (3^:^^(1-0''^*) and k := k - ^p. 

4k 

Proof. Under the stoppcd-sharc-price measure (|2.18p the Heston dynamics arc given by 

dA:„ = (- 114 + Vulu<t) du + VKdWu , Xo e M, 

(5.27) dVu = {kO - nVa + p^VuluKt) du + i^dZu, Vq ^ v > 0, 
d{W,Z)^ ^pdu. 

We now compute 

(5.28) E (^c"(^'+--^')) = E (e (c"(^'+--^')|J-t)) = E (^c'^M+Biu.rW,^^ 

(5.29) = e^("'^'E (c^(".^)^' j = e^("'^'A^(B(u, r)). 



30 



ANTOINE JACQUIER AND PATRICK ROOME 



where 

(5.30) AY{u) = exp I ^ j (l - 2/3tu) for all " < 



with q := 4^6*/^^, A := ue~'^*/A- Over + t] the Heston dynamics in (|5.27p remain the same as under the 
risk- neutral measure and so we can apply the standard spot Heston mgf in (j5.28p . On [0, t] we use the modified 
chi-squared mgf in (|5.30l) to obtain (|5.29p . corresponding to the modified Heston process in (j5.27p on [0,t]. □ 

5.2.2. Proofs of Section [S.1.2[ In this section we shall use the standing assumption that k > p^. Let e = r^^ 
and consider the Heston process (r-i^l*') with (Xt)t defined in and defined in (|2T1|) . 

Specifically defined in (|2.ip is then given by A.(u) = T-^E fe"-^'") , and for ease of notation we set 



(5.31) A(.*)(w) = Ae(w) forahwePe- 

The following lemma [501 page 13] recalls some elementary facts about the function V in (|3.12l) . which will be 
used throughout the section. We then proceed with two technical results needed in the proof of ProDOsition l5.13l 

Lemma 5.9. The function V defined in p.l2p is infinitely differentiahle, strictly convex, essentially smooth on 
the open interval {u^,u^) with u± defined in p.Sp and 



k9 {2k — + prf) 

u\u\ ■ 2^2 (1 _ ^2) 



V(u^) := hm V(u) = ^ < oo, 



V (u+) hm V{u) = La// 2\ < 

u7'«+ 24^ (1 - p'^) 



with u_ < 0, > 1 and V{0) = V{1) = 0. Furthermore, it has a unique minimum at 2^{i-p^) ^ 

Lemma 5.10. Let p± be defined as in p. Ill) . /3t in p.l3p . and recall the standing assumption p < k/^. 
Assume further that t > and define the functions g+ and by 



9± (p) := (2^ - pO ± p^e{l-p') + {2n-pO' ^ ^\ " ^ 

Pt 

(i) The inequality —1 < < always holds and if k/^ > then 1/2 < < 1; 
(ii) the inequality g+{p) > holds if and only if p £ {p+, 1) and k/^ > p+; 
(Hi) the inequality g-{p) > holds if and only if p (z {—l,p-). 

Remark 5.11. From the proof, we have the equality p_ — —1 if and only if t = 0. Also if k > then 
/9+ = 1 if and only if t = 0. 

Proof. We first prove Lemma [5.10r i) and consider the inequality — 1 < p_ < 0. Using the definition of p_ 
in p. lip this is equivalent to proving that 

^ - + f < - Vl6.2e2^^ ^eH- e-)- < ^ (1 - e-) = ^^j ' ■ 
e''* + 1 V ^ ^ ^ s V y 1 + e"^* 

The upper bound then follows trivially. Also we can write 



16K2e2«t+^2(l-e«*)^ = 



2_ / (■? - (8k + Oc^"*)^ 16Ke2'^* (e''* - 1) (k + C + ^e''* + 3Ke«*) 



(e«* + 1)^ (e-^* + 1)^ 
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and the lower bound follows. We now prove that p+ > 1/2. From (|3.1ip this is equivalent to 
Tedious rearrangements show that the left-hand side ean be written as 



16^262'^* +C2(l-e«*)^ = 



' (4^ +{k- 40e2«*)^ Ke2'^* (8^ (e2«* - 1) + n (512e''* + 255e2''* + 256)) 
16 (e«* + if 16 (e'^* + 1)^ 



and the result follows. We now prove that the upper bound p-^- < I holds if k/S, > p+. Assume that n/S, > p+. 
Since g^^^^t^^ is always strictly smaller than 1/2, we immediately obtain the inequality 

Using the definition of in p. lip the statement p+ < 1 is equivalent to 

e + (8k - C)e2'^* 



lg^2e2«:t_^^2(l_eKt)2 < 



which can be written as 



/ + (8k - g)e2«*)2 16Ke2«* (c"^* - 1) (k - ^ je'^* + 1) + Skc'^*) C + (8k - g)e2'^* 
V (e''* + 1)^ (e''* + 1)^ ^ e«* + l ' 

This statement is true if k — ^ (c'^* + 1) + Skc"* > and if the rhs is positive. The former inequality is 
precisely (|5.32p . and immediately implies the positivity condition. Therefore p+ < 1. 

We now prove Lemma l5.10r ii). The equation g+{p) = implies (by squaring and rearranging the terms): 

(5.34) 4K(p2 _ 1) (4Ke2'=*p^ + ^(1 - c^'^')^ - k(1 + 20*^* + e2'=*)) = 0. 

The roots of this equation are ±1 and p± defined in (|3.1ip . Clearly some of these solutions are extraneous. The 
two possible positive roots are 1} and the two possible negative ones are —1}. Clearly g+(— 1) = 0. 
Straightforward computations show that g'j^{—l) < and g'^{0) > 0. Since is continuous on (—1,0) with 
g+{0) < 0, it cannot have a single root in this interval, and p_ £ ( — 1, 0) (by Lemma lS.lOr i)) is hence not a valid 
root. Consider now p £ (0,k/^). From Lemma [5.101^ 1) the only possible roots are p £ {l,p+}. If p = 1 then 
clearly k/^ > 1 by the standing assumption. But .g+(l) = 2k — ^ + ^/(2K^^^y^ is null if and only if k/^ < 1/2, 
which is a contradiction. Therefore the only possible positive solution is p+. Now, on (0, k/^) we have 

9'4p) = -c - , , ^, + V4k2 - 4Kep +e + 

^J^K? - \k^p + ^2 1 - e 

>-e- ... '"5:^ .. + V4K2-4KaK/e)+g2+ 



8k/9 



2kp > 0. 



1-e- 

In summary, g+ is strictly increasing on (0, k/^) with a unique zero at p+ satisfying k > jO+^. On the inter- 
val (— 1,k/^), g+{p) > if and only if p e (p+jl) and p+(, < k. The proof of (iii) is analogous to the proof 
of (ii) and we omit it for brevity. 

□ 

Lemma 5.12. Let p± and u^. be as in (j3.1ip and (j3.8p and t > 0. Then u*^ > ^ if P l£ P-- o-nd u*_ < Q if 
P > P+- 
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Proof. From p.9p write v = \J z{p), where 

The two numbers u*_ and in p.Sp are well-defined in M if and only if z{p) > and t > 0. The two roots 
of this quadratic polynomial are given by x± '■= 37: [e"""* (^(c"* ~ 1) i inc'^*/^^] . We now claim that p_ < x- 
and > From the expression of p_ given in p. lip , the inequality p_ < x- can be rearranged as 



Noting the the square root term is equal to 



^ - 2^e''* + ^e2«* 



8/te 2 



4eKt (e^t _ 1)2 / ^ ^ 2^0"^ ) ( C - 2Ce'^* + ^e- 



(gKt _^ 1)^ (e«t _^ 1)^ 

proves the claim. Analogous manipulations imply the other claim > x+j 8.nd hence z{p) is a well-defined 
real number for p G [— U 1]. 

The claim u*_ < is equivalent to -^^2 _ 2e''* (8^2 - 4k^p + ^2) _|_e2Kt(-^ _ 4^5^)2 < ^ (1 _ e"*) 4- Anpe*^^, 
which holds as soon as ^ (1 — e"*) -|-4Kpe'** >0, orp> ^(1 — e^***). Therefore the claim follows for any p > p+ 
if and only if p+ > ^ (1 — e^'^*). This inequality simplifies to 



V^2 + 16«;2e2«* - 2C2e«t + ^2g2«t > 



C(e''*- 1)' 



1 ' 



which can be written as 



4e''* (^4K2e''* (e''* -1- 1) 


2+^2 (e^t_ 1)2^ 


1 ,^2(eKt_i)4 


(e«* 4- 


If 


(e«* 1)^ 



> 



1 ' 



which is always true. 

Now straightforward manipulations show that the inequality u*^ > 1 is equivalent to 



which is true if p < 



(5.35) 



(e«* Anpe'^^Y - 16kc''* {k + ^p (e«' - 1)) > C (e''* - l) + 4Kpe'"*, 

^(1 _g-«t) 



£, (e«* - 1) 



or p < 



(5.36) 



4k 
P-<- 

or p_ < — 



And of course the claim > 1 if p < p_ ) shall hold if 



^ (e«* - 1) 
C(l-e-''*) 



An 



Consider (|5.35p . This inequality, which can be re- written as 



15^2e3Kt ^^2 (-gKt _ 1)2 (gKt + 1) - 4k2c''*^ 


/^2(1 _ g,ct)(l _ g2«t) ^ 8^2g2«t\ 2 


^2 (g2,.t _ 1)2 ^ 


C(e''* + l)(l-e'"*) J 



< 



g2(l_eKt)(i_g2.:t)_^g^2g 

C(e«* + 1)(1 - e«*) 



,2„2Kt 



holds if (e-^* - 1) (e-^* + 1) - 4K2e'"* > 0, or 



> TT- Quick manipulations turn ()5.36p into 



4gKt |^4K2gKt (gKi _^ 1)2 _ ^2 (gKt _ 1)2 (2e«^* + 1)^ 


^2 (2e«* - 3e2'^t + 1)2 


(e«* + 1)^ 


(e«* 1)^ 



< 
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Again this trivially holds if 4K^e'** (e*** + 1)^ - (e''* - 1)^ (2e'^* + 1) > 0, which is in turn equivalent to |^ > 

-^^ — ^^^t\ly^ — ■ omce the inequality -^^ — 4(c'^t+i)^ — ^ 4 ' clearly true, it follows that for any 

valid choice of parameters either (|5.35p or (|5.36p (or both) hold, and the claim follows. □ 

Proposition 5.13. Let e ~ t^^ and consider the large-maturity Heston forward process (t-^x'tA withXt 

\ /t>0 

defined in p.f p and X^*^ defined in (|2.f fp . Then I?o.o = ^oo ^.n-d {0, f } C Vq q with 'Doc, Dq q and 2?q q defined 
in p.lOp and in Assumvtion \2. 1\ 

Proof. We write 

E (e"(^'+--^')) = E [e (e"(^'+--^')|J-t)] = E (^e^K-)+s(n,r)v,^ ^ e^^^^^^E (e^("^-)^*) . 
For any fixed t >Q we require that 

(5.37) E (^e"(^*+--^')|J"t) < oo for ah t > 0. 

Andersen and Piterbarg |31 Proposition 3.f] proved that if the following conditions are satisfied 

(5.38) K > p£_u, 

(5.39) (k- p^uf + u{l~u)^'^ >0, 

then the explosion time is infinite and (|5.37p is satisfied. In [50] the authors proved that these conditions are 
equivalent to k > p£, and u G with w_ < and u+ > I {u± defined in p.8p ). Further we require that 

(5.40) E (c^f"'^)^*) < oo, for aU r > 0. 



Now denote 

Vv ■■= {ueR-.E (^c^i^'-^)^t^ < oo, for ah r > 0}. 

Then the limiting forward mgf domain is given by I?oo = ['u_,it+] H Vy and k > p^. The condition (|5.40p is 
equivalent to B {u,t) < l/(2/3j) for all r > 0, where 

p. N £-2/ c J/ 1 - exp (-d (m) r) 

1 — 7 (u) exp [—d (u) T) 

and where d and 7 are given in p.lSp . Now a simple calculation gives B (0, r) =0 and _B (1, t) = for all t > 0. 
Furthermore for u £ (0,1), and given Conditions (|5.38p and (|5.39p . we have d{u) > k — p^u and 7 (u) < 0. 
This implies that B {u,t) is strictly negative for u e (0, 1) and t > 0. In particular [0, 1] C I?oo (recall that the 
process is a martingale). For a fixed u G M we calculate 

dB {u, t) _ 2u{u - l)d(u)2e''(")^ 

dr (^ _ i^ed(u)T _^ l^^d{u)r - 1) - d{u) (e''(")^ + 1))^ ' 

so that for any u ^ [0, 1], B {u, •) is strictly increasing. Therefore the limiting domain is given by 

2?oo = |w e M : lim S(u,t) < 1 H 

L r^oo 2^t J 

with [0, 1] C I?oo- Wc first concentrate on the first condition and check at the end that our solution always 
contains [0, 1]. We have limT-^oo B [u, r) = [k — p^u — d (u)) = V{u)/{k9), with V defined in p.l2p . So the 
condition is equivalent to k — p£^u — d{u) < 2k/(1 — e"*^*). If p < (p > 0) and u < {u> 0) then 

2k 

K — p£,u — d{u) < K — p£,u < K < —7, 
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and the condition is always satisfied. So if p = the domain is given by . If p < (p > 0), the 

domains contain [m_,0] ([0,u+]). Now suppose that p < and u > 0. The condition above is equivalent to 
V{u) < K9/{2/3t). From Lemma [5.91 on (0,m+], the function V attains its maximum at u+ with 

kO {2k -pi,- pij) 

Using the properties in Lemma l5.9| there exists m^J. G (1, u^) solving the equation 

(541) mi^± 



if and only if (2k — p^) — p\J (1 — _|_ — p^)^ > (l — p^) //3f. This condition has been shown in 
Lemma r5.10l to be equivalent to — 1 < p < p_ and t > Q. The solution to (|5.4ip has two roots u*_ and u*^ defined 
in p.Sp . and the correct solution in this case is u*^ since from Lemma 15.121 we know that u*^ > 1 and ul < 0. 
So if p_ < p < then the effective limiting domain is . If — 1 < p < p_ and t > Q then the effective 

limiting domain is given by [u_ , it!j_) . Following a similar procedure we can show for < p < min (k/^, p+) the 
effective limiting domain is given by . If p_|_ < p < min(K/^, 1), t > and ti > p+^ then the effective 

limiting domain is given by (itl,M+] , with <u*_ <Q. □ 

Lemma 5.14. The following expansion holds for the forward mgf K^t^ defined in ()5.3ip ; 

A(.*'(u) = V{u) + ^^-^ {l + O (^c"'^("'^)) , for all ueV^o, as t tends to infinity, 

where the functions V , H, d and the interval Voo are defined in p.l2p . (j3.13p and (j3.10|) . 

Remark 5.15. For any u e I?oo, d{u) > 0. Indeed I?oo ^ [w-jU+J by Proposition 15.131 Furthermore since 
K > p^, u S u+] implies (j5.39p which in turn implies d{u) > 0. 

Remark 5.16. We note in particular that the exponential decay in the remainder implies that for all u € Doo, 
Ar '(u) = V{u) + H{u)/t + O (I/t'^) for r large enough, which is used in the proof of Proposition 13.81 

Proof of Lemma \5.14\ From the definition of in (|5.3ip and the Heston forward mgf given in (|5.16p we 
immediately obtain the following asymptotics as t tends to infinity: 

A (., .) . rViu) - ^ log (^) + O (e-^").) , Biu,r)^^ + (e-(«)^) , 
where A and B arc defined in (|5.17p . V in p. 121) and d and 7 in (|3.13p . In particular this implies that 

l-2/3tB{u,T) eK-2/3tV{u) \ J' 
log (1 - 2P,B {u, r)) = log (1 - + o (e-'^(-)-) , 

which are well-defined for all u £ I?oo- We therefore obtain 

^, , Viu) 2^9^ / 2l3tV{u)\ 2kB^ 

= -TT^TTTTTX^'^ " log U ITT^ - log 



and the lemma follows from straightforward simplifications. □ 
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5.3. Proofs of Section 13. 2i For ease of notation here, the subscript i for a function or a variable shall refer 
to the i-th Hcston model, for which we can readily use the results of Section [5TTJ For instance the function A; 
shall refer to the function A evaluated using the parameters in the i-th Hcston model. 

Lemma 5.17. In the n- dimensional multivariate Heston model (|3.15p the forward mgf defined in (|5.15p reads 

n n 

A{u) = logE (c"-^-") = ^ Ai(w), for all u e f] 2?a,, 

1=1 i=l 

where Ai and 'D\. are defined in (j5.16l) . 

Proof. Conditioning on the filtration (J^t)t>o and using the tower property, we obtain 

n 

A{u) = logE (e (c"^''>t)) - ^ A,(7.,r) +logE (c^r^i ^'(".-)^/'') , 

i=l 

using the mgf of the standard Hcston model. By independence of the variance processes, we obtain 

n n n 

A(u) = ^A,(w,T)+logE(c^"=i^'("^^'^''") =^A,(w), forallwe p|2?A,. 

□ 



i—1 i—1 2—1 



5.4. Proofs of Section 13.31 Consider the following functions: 
(5.42) 

Aiu, r) :^ .) + ^"'^'(xW - ^(~)) J,,„, _ , _ ^ , ,„-.,*5iWiMp^*-"' 

A I \ 1/ / \ 1/ \\ li /I - 7 (w) exp (-2d(w) t) 
^1 (w, r) := - - d (w)) T - - log ' 



2 ''^^ ' " " 2 °\ 1 - 7(u) 

AKe x{u)~d{u) (l-cxp(^dHT))' 2(xH-d(u)) 1 - exp (-2d (u) r) 

Si(w,r):=— 5 IT^^ -; r^rrr^TT' ^2(m,t) — 



d(u) 1 — 7(u) exp (— 2(i(ii)T) 1 — 7 (m) exp (— 2d (w) r) 

and 



(5.43) M(r,p, log (l-2rg)j, |- (l - e^^Kt) ^ /x* := V^^e-^' + (l - e"'^*) , 

Note that although in some cases we use the same variables and functions as in the Heston analysis they 
may have a different definition in this section. In our analysis we require the following lemma which is a direct 
consequence of [U Equation 29.6]. 

Lemma 5.18. If Z ^ J\f{0, 1) and p and q are two constants, then, for AI defined in (|5.43p . 

logE (^e"(P^+''^')) = M{u,p,q), whenever uq < 1/2. 
Lemma 5.19. In the Schohel-Zhu model (j3.18p the forward mgf defined in (|5.15p is given by 

A{u) = A{u, t) + Bi{u, T)nt + B2{u, t)^4 + m(i, Bi{u, T).JJt + 2B2(u, T)^/Jtl^t, B2{u, r)/?*) , 
for all u e I?A, where all the functions and variables are defined in (|5.42p and (|5.43p . 
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Proof. Conditioning on the filtration {J-'t)t>a and using the tower property we find 

A(w)=logE E (^e"-^*" I J"t) ^ A{u,t) +logE[cxp{Bi{u,T)at + B2{u,t)(t'})] , 

where we have used the standard Schobel-Zhu mgf from [23 and the functions defined in (|5.42p . Since at is 
Gaussian with mean fit and variance /3t (given in (|5.43p ). we obtain 

A(ii) A(M,r) + logE (e^i("^^''"'+^^("^^''"') 

with Z ^ Af{0, 1), and the lemma follows directly from Lemma [5. 181 □ 

5.5. Proofs of Section 13.41 Let (j> be the Levy exponent of the Levy process Y. If v follows (|3.2ip . by a 
straightforward application of the tower property for expectations, the forward mgf defined in (|5.15p is given by 

(5.44) A(^.) = A(0(u), r) + —^IMlI) „e~'=* - ^ log (1 ~ 2/3ti3(0(u), r)) , for all u e Pa, 
where 

k6 f , ,, , / 1 — 7(m) exp (— d(?i)T) 

(5.45) A{u,T):= — [{K-d{u))T-2\og' 

B{u,t) 



e V V i-7(") 

K — d{u) 1 — exp (— d(it)T) 
^2 1 — 7(u) exp (— (i(u)r) ' 



(5.46) d{u):={n'-2uef\ ^("^ ^T^Sl '"""^ ^* ^= £ 

Similarly if {vt)t>o follows (|3.22p the forward mgf defined in (j5.15p is given by 

^A* , f B{<t>{u),T 

^et\B{c^{u),T)-a) 



(5.47) A{u) ^ A(0(M),r) +g(,^(^),T)^;e"^* +dlog ( ^^y^^ilJ for aU u G Pa, 



where 



Ad 



u 



(5.48) A(u,t)-~— ur + alog 1 r (l - e"^"") and B(u,t):=-(1 



u 



A 



e 



-Xt\ 



Proof of Proposition \3.19[ We show that Proposition 12.141 is applicable given the assumptions of Proposi- 
tion [3?lll Consider case (i). The expansion for A^'' defined in (|5.3ip is straightforward and analogous to 
Lemma 15.141 In particular we establish that 



AW [u) ^ V{u) + ^(^l + (e-'^f")-) ) , 



for all u G Pqo , as r tends to infinity, 



where the functions V , H, d and the domain I?oo are defined in p.23p . (|5.46p and (|3.25l) . Since 4> is essentially 
smooth and strictly convex on and I?oo ^ V^, then the limiting mgf Ao,o = V is essentially smooth and 
strictly convex on Vao- Also Ar is infinitely differcntiable on Poo since is of class C°° on Poo- Since (/>(!) = 
we have that V'(l) = and {0, 1} C P^. It remains to be checked that the limiting domain is in fact given 
by Poo- We first note that that by conditioning on {Vu)t<u<t+T and using the independence of the time-change 
and the Levy process we have E (c"^^''+^~^''^) ~ E ^c''^^"'-'*^ ^^ds^ ^^^^ ^^^-^y ^ ^^ic limiting domain must 
satisfy (/)(u) < oo. Using [Ml page 476] and the tower property we compute 

(5.49) E (c"(^^*+^-^*)) = E [e (e^(«)/r^''='i^|J-t)] = E (e^('^(«).-)+B(<^(«).r)., j ^ ^Ai4>iu),r)^ (^^BW^),r)v,\^ ^ 
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with A and B defined in (|5.45p . Further from (|5.18p we have 

logE (e™* ) = - — log (1 - 2(3,u) , for all < — . 

1 - 2ptu C 2ft 

Following a similar argument to the proof of Proposition 15.131 we can show that for any t > 0, B{(j){u),T) < 
l/(2/3t) is always satisfied for each t > 0. This follows from the independence of the Levy process Y and the time- 
change. We also require that for any t > 0, E ^e-^^* ^ "sds0(u) | jr ~j ^ every r > 0. Here we use [5] [Corollary 
3.3] with zero correlation to find that we require 4>{u) < k^/(2^^). It follows that = {u ■ </'(u) < '^^/(2C^)}. 

Regarding case (ii), arguments analogous to case (i) hold and we focus on showing that the limiting domain 
is Vrx,- Using [T4j page 488] Equality (|5.49p also holds with A and B defined in (|5.48p . Since we require that 
for any t > 0, E (^e-^^t*^ ^"'^'^'^^^'^{Tt^ < oo, for every r > we have < aX. Using [T31 page 482] we also have 

f 11 - — rvp^^ \ 

log E (e""' ) = e-^* + log — , for aU u < a. 

\[u — a)e^' / 

But it is straightforward to show that < aX implies B{(j){u),T) < a for any r > and it follows that 

T^oo — {u ■ < aX}. Case (iii) is straightforward and omitted. □ 
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